m 


la.  REPORT  SEC 


2a.  SECURITY  CLASSIFICATION  AUTH 


2b  DECLASSIFICATION  /  DOWNGRAD 


-  1  gBg 


6a.  NAME  OF  PERFORMING  ORGANIZATION 


6c  ADDRESS  (City.  State,  and  ZIP  Code) 
Pittsburgh,  PA  15260 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


8c  ADDRESS  (City,  State,  and  ZIP  Code) 
AFOSR/NM 
Bldg  410 

Bolling  APB  DC  2 


1 1 .  TITLE  (Indude  Security  Classification ) 

Multival^^ate  Analysis  and  its  Application 


12.  PERSONAL  AUTHOR(S) 

P.R.  Krishnaiam  &  C.R.  Rao 


13a.  TYPE  OF  REPORT 

Final 


16.  SUPPLEMENTARY  NOTATION 


1b  RESTRICTIVE  MARKINGS 


3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release  ; 
distribution  unlimited. 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

AFOSR.TR.  8  7-1 


6b  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 
(If  applicable) 

AFOSR/NM 


7b  ADDRESS  (City,  State,  and  ZIP  Code) 
AFOSR/NM 
Bldg  410 
Bolling  AFB  DC 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 

61102F 


PROJECT 

TASK 

NO 

NO 

2304 

A5 

WORK  UNIT 
ACCESSION  NO 


13b  TIME  COVERED 

FROid0/l/84  to  9/30/87 


14.  DATE  OF  REPORT  (Year,  Month,  Day)  IIS.  PAGE  COUNT 
Sept.  1987  I 


COSATI  CODES 


SUB-GROUP 


18  SUBJECT  TERMS  ( Continue  on  reverse  if  necessity  end  identify  by  block  number) 


19  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

This  contract  supported  the  Center  For  Multivariate  Analysis  which  included  support  for 
two  core  researchers  and  several  visiting  researchers.  A  total  of  160  technical  reports 
were  issued;  of  these  48  have  been  published,  38  have  been  accepted  for  publication,  and 
74  have  been  submitted  for  publication.  Some  important  contribution  were  made  to  several 
areas  of  multivariate  analysis,  with  applications  in  manufacturing  technology,  signal 
processing,  automation,  expert  systems,  pattern  recognition,  and  machine  intelligence. 


20  DISTRIBUTION/ AVAILABILITY  OF  ABSTRACT 

□  UNCLASSIFIED/UNLIMITED  □  SAME  AS  RRt  □  3T'C  USERS 

21  ABSTRACT  SECURITY  CLASSIFICATION 

22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Maj .  Woodruff 

22b  TELEPHONE  (Include  Area  Code) 

(202)767-5027 

22c.  OFFICE  SYMBOL 

NM 

00  Form  1473.  JUN  86 


Previous  editions  are  obsolete 


5ECURITY  CLASSIFICATION  QF  This  PAC-E 


r  V  ..v.  s; V. j, v V.v.v  . 


•a  .3  4  AT  <* 

/.V  V  V  ’ 


vssv, 


AF08H.TR 


8 7 • 1 967 


mwwwwwwwwww 


MULTIVARIATE  ANALYSIS  AND  ITS  APPLICATIONS 
Final  Report  Under  the  Contract  F49260-85-0008 

Sponsored  bv 

the  Air  Force  Office  of  Scientific  Research 
October  1984  -  September  1987 


P.R.  Krishnaiah  and  C.R.  Rao 
Principal  Investigators 


i 


This  is  the  final  report  on  the  research  work  done  during  the  3-year  period  October 
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Part  1.  SOME  HIGH  LIGHTS  OF  THE  RESEARCH  WORK  DONE 

During  the  3-year  period,  some  important  contributions  were  made  to  several  areas 
of  multivariate  analysis.  They  will  be  specially  useful  in  solving  some  problems  of 
interest  to  the  Department  of  Defense.  For  example,  they  have  applications  in 
manufacturing  technology,  signal  processing,  automation,  expert  systems,  pattern 
recognition,  machine  intelligence  etc..  A  brief  description  of  some  of  the  significant 
contributions  is  given  below. 

1.1  SIGNAL  PROCESSING 


When  a  message  is  received,  it  may  be  pure  noise  or  it  may  be  a  specified  or  an 
unknown  signal  mixed  with  noise.  The  problems  that  arise  are  detection  and  extraction 
of  signals  from  a  rece;.ved  message.  Satisfactory  solutions  are  obtained  in  the  following 
situations. 

1.  A  test  is  developed  to  decide  whether  a  received  message  contains  a 
specified  signal  (including  the  case  of  no  signal)  when  the  noise  distribution 
is  Gaussian,  elliptically  symmetric  and  so  on.  The  test  is  based  on  an 
estimated  covariance  matrix  for  noise,  and  its  power  depends  on  the 
degrees  of  freedom  of  the  estimated  covariance  matrix.  Some  guide  lines 
are  provided  for  the  choice  of  degrees  of  freedom  to  ensure  a  desired  level 
for  the  power  for  a  specified  signal. 

2.  A  model  for  a  received  message  which  has  been  extensively  studied  is  of 
the  form  x(t)  =  A  s(t)  +  n(t)  where  n(t)  is  pxl  noise  vector,  s(t)  is  qxl  signal 
vector  and  A  =  [ A(<J>j  ),■••.  A(pq)]  is  pxq  matrix  whose  columns  are 

functions  of  unknown  parameters.  Under  the  above  m  del,  methods  are 
developed  for  estimating  the  number  of  signals  and  also  the  signal 
parameters.  The  problem  of  estimation  of  the  number  of  signals  is  treated 
within  the  framework  of  model  selection  procedures  using  some  new 
information  theoretic  criteria.  The  problems  are  solved  under  white  and 
colored  noise  distributions  for  the  errors.  The  new'  methods  are  more 
general  and  flexible  than  those  that  are  currently  in  use. 

3.  Another  model  in  the  area  of  signal  processing  is  of  the  form 

x(t)=u lexp(io)lr)+  +d^exp(i<Si(Jt)+n(t),  r=l _ V 

involving  multiple  sinusoids.  The  parameters  ar  ov  and  q.  the  number  of 
signals,  are  unknown.  The  usual  non-linear  least  squares  theory  and  tests 
involve  heavy  computations.  A  new  method  of  estimation  called  EVLP 
(equivariation  linear  prediction)  has  been  proposed  to  estimate  the  number 
of  signals,  amplitudes  and  frequencies  of  signals.  The  estimates  obtained 
by  EVLP  have  nice  properties,  are  easy  to  compute  and  are  more  efficient 
than  some  of  the  available  methods. 

4.  In  the  same  context  as  in  (2)  and  (3)  above,  and  in  a  more  general  situation 
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of  superimposed  exponential  signals,  some  non-parametric  procedures  have 
been  proposed  to  estimate  the  number  of  signals.  They  are  based  on 
general  information  theoretic  considerations. 

5.  Alternative  criteria  based  on  Loo  (leave  one  out)  method  have  also  been 
developed  for  the  estimation  of  number  of  signals.  This  method  is  easy  to 
implement  and  is  useful  in  small  samples. 

1.2  DENSITY  ESTIMATION 

The  discriminant  function  in  a  classification  problem  is  based  on  the  ratio  of 
probability  densities  in  two  alternative  populations.  In  practice  one  has  to  estimate  these 
densities  from  observed  samples  from  each  population.  In  a  parametric  situation,  the 
estimation  of  the  density  function  is  done  through  the  estimation  of  the  unknown 
parameters.  In  a  non-parametric  situation  there  are  several  methods  available  in  the 
literature  on  density  estimation  in  the  univariate  case.  Special  methods  have  been 
developed  in  the  multivariate  cases  when  we  have  binary  data  and  multidimensional 
directional  data.  The  latter  is  useful  when  we  want  to  discriminate  using  only  the 
information  on  shape  (ratios  of  measurements). 

1 .3  CHANGE  POINT  PROBLEMS 

Problems  of  detection  of  change  points  arise  in  quality  control  programs.  For 
example,  in  continuous  production  it  is  important  to  find  out  when  a  deterioration  in  the 
quality  of  a  product  takes  place.  They  also  arise  in  problems  of  edge  detection.  Studies 
have  been  made  on  estimating  the  number  and  location  of  change  points.  The  strong 
consistency  of  the  proposed  procedure  is  established. 

1 .4  ANALYSIS  OF  CONTINGENCY  TABLES 

The  well  known  likelihood  ratio  and  chi-square  tests  for  independence  of  attributes 
in  contingency  tables  are  not  efficient  when  the  categories  of  each  attribute  are  known  to 
be  ordinal  in  namre  and  the  alternative  hypothesis  is  specified  to  be  a  particular  type  of 
dependence. 

Some  extreme  point  methods  have  been  introduced  in  analyzing  problems  of  this 
type.  These  methods  are  useful  in  obtaining  precise  expressions  for  the  size  and  power  of 
alternative  tests.  As  a  consequence,  it  is  much  easier  to  select  a  test  depending  on  a 
particular  problem. 
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1.5  GROWTH  CURVE  MODELS 


The  problem  of  predicting  a  future  measurement  on  an  individual  given  the  past 
measurements  has  been  studied  using  a  growth  curve  model  for  the  measurements.  The 
efficiencies  of  different  methods  of  prediction,  such  as  direct  regression,  inverse 
regression,  extrapolating  a  growth  curve  fitted  to  an  individual's  past  measurements, 
empirical  Bayes,  factor  analytic  regression  etc.,  are  examined  by  cross-validation  or 
leave-one-out  technique.  A  number  of  data  reduction  methods  are  developed  and 
analyses  based  on  them  are  discussed.  A  new  method  known  as  calibrated  prediction  is 
developed. 

1.6  MULTIVARIATE  CALIBRATION  PROBLEMS 

The  problem  of  predicting  response  variables  using  explanatory  variables  is  well 
known  and  this  is  accomplished  by  computing  the  regression  of  the  former  on  the  latter. 
But  situations  arise  where  it  is  of  interest  to  determine  the  values  of  the  explanatory 
variables  knowing  the  values  of  the  response  variables.  Such  a  setup  is  called  calibration. 
The  multivariate  calibration  problem  can  be  described  theoretically  through  the  equation 

X-SH+Bx  +  e 

where  y  is  the  vector  of  response  variables,  x  is  the  vector  of  explanatory  variables,  B  is 
the  matrix  of  regression  coefficients  and  a  is  constant  vector.  Methods  have  been 
developed  under  the  condition  that  e  ~  Np(0.Z),  i.e..  e  is  distributed  as  p-variate  normal 
with  mean  zero  and  covariance  matrix  Z.  for  estimating  x  given  y  and  for  selecting  the 
best  subset  of  y  for  predicting  x. 

1.7  TIME  SERIES  ANALYSIS 

New  methods  have  been  developed  to  estimate  the  order  (.moving  average  and 
autoregressive  terms  >  of  an  ARMA  process  based  on  time  series  data.  The  building  of  an 
ARMA  model  applicable  to  a  given  time  series  data  is  viewed  as  a  search  for  a 
parsimoniously  parameterized  linear  filter  that  transforms  the  observed  time  series  into 
white  noise,  i.e..  a  icm  e  of  uncorrelated  random  variables.  Model  detection  criteria 
are  used  to  test  tor  aite  noise.  which  automatically  assigns  a  significance  level  for  the 
test,  which  depend'  on  the  sample  size  and  the  model  detection  criterion  used  such  as  the 
AIC.  BIC  and  so  on  Extensive  simulation  results  have  shown  the  superiority  of  the 


Part  2:  LIST  OF  TECHNICAL  REPORTS  AND  ABSTRACTS 

All  the  Technical  Repons  were  written  with  partial  or  complete  support  under 
contract  F49620-85-0008  with  the  Air  Force  Office  of  Scientific  Research  during  the 

period  October  1,  1984  -  September  30,  1987 

1.  Bai,  Z.D.,  Limiting  Properties  of  Large  System  of  Random  Linear 
Equations.  Technical  Report  #84-41, Center  for  Multivariate  Analysis, 
October  1984. 

S.  Geman  and  Chi  R.  Hwang  proposed  a  kind  of  algebraic  system  of 
equations  and  proved  the  law  of  large  numbers  for  its  solution.  In  this 
paper,  the  conditions  to  ensure  these  results  are  significantly  weakened  for 
the  law  of  large  numbers.  Also,  the  central  limit  theorem  is  shown.  For 
both  the  law  of  large  numbers  and  the  central  limit  theorem,  the  only 
needed  assumption  is  that  the  random  variables  have  finite  second  moment. 

2.  Bai,  Z.  D.,  Yin,  Y.  Q.,  and  Krishnaiah,  P.  R.  On  Limiting  Empirical 
Distribution  Function  of  the  Eigenvalues  of  a  Multivariate  F  Matrix. 
Technical  Report  #84-42,  Center  for  Multivariate  Analysis,  October  1984. 

In  this  paper,  the  authors  derived  an  explicit  expression  for  the  limit  of  the 
empirical  distribution  function  of  a  central  multivariate  F  matrix  when  the 
number  of  variables  and  degrees  of  freedom  tend  to  infinity  in  certain 
fashion.  This  distribution  is  useful  in  deriving  the  limiting  distributions  of 
certain  test  statistics  which  arise  in  multivariate  analysis  of  variance, 
canonical  correlation  analysis  and  tests  for  the  equality  of  two  covariance 
matrices. 

3.  Bai.  Z.  D.  The  Strong  Consistency  of  Error  Probability  Estimates  in  NN 
Discrimination.  Technical  Report  #84-43.  Center  for  Multivariate 
Analysis.  October  1984. 

A 

Since  1971.  T.  J.  Wagner  introduced  the  error  probability  estimate  Rn  in 
NN  discrimination  rule  and  proved  the  weak  convergency  of  this  estimate. 

In  this  paper,  we  show  the  exponential  convergence  rate  of  this  estimate. 
Hence,  the  strong  consistency  of  Rn  holds  automatically. 

4.  Bai,  Z.  D..  Krishnaiah.  P.  R.,  and  Yin.  Y.  Q.  On  Limit  of  the  Largest 
Eigenvalue  of  the  Large  Dimensional  Sample  Covariance  Matrix. 
Technical  Report  S84-44.  Center  for  Multivariate  Analysis,  October.  1984. 

In  this  paper,  the  authors  showed  that  the  largest  eigenvalue  of  the  sample 
covariance  matrix  tends  to  a  limit  under  certain  conditions  when  both  the 
number  of  \  enable-;  end  the  sample  si^e  tend  to  infinity.  The  above  result 
is  proved  unde:  tire  mild  restriction  that  the  fourth  moment  of  the  elements 
of  the  sample  of  squares  and  cross  products  iSP'i  matrix  exist. 

5.  Fang.  Zhaob-.n  \  Note  about  the  Strong  Convergence  of  the 
N’onparatr.e:-;  !  notation  of  a  Regression  Function.  Technical  Report 
#84-45.  Center  t . >r  Multivariate  Analysis.  September  1984. 

Consider  the  regression  model  Y,*  =  g(x.*)  +  e  *  .  i=1.2.  ,n.  x,*'s  are 
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unordered  design  variables,  g  unknown  function  defined  on  [0,1].  {e^} 
i.i.d  r.v  with  mean  0  and  finite  moment  of  order  p  >  1.  The  asymptotic 
behavior  of  estimator  ^  are  studied. 

6.  Bai,  Z.  D.,  and  Yin,  Y.  Q.  Limiting  Behavior  of  the  Norm  of  Products  of 
Random  Matrices  and  Two  Problems  of  Geman-Hwang.  Technical  Report 
#84-46,  Center  for  Multivariate  Analysis,  November  1984. 

Let  Wn  =  (w-)  i.j=l,2,  ■  •  ,n  be  a  matrix  of  iid  entries  with  mean  zero, 

variance  c~  and  finite  fourth  moment.  The  first  result  of  this  paper  is  that 
W 

"tim  II  — _  II  <  (l+k)Gk,  a.s.  where  II  •  II  denotes  the  operator  norm  in  L-,. 


As  a  consequence  of  this  result  we  have  that  tim  SR{— r)<G,  a.s.  ,  where 

"  Sn 

SR  stands  for  spectral  radius.  Applying  the  first  result,  we  solved  two  open 
problems  of  Geman  and  Hwang  concerning  the  asymptotic  normality  of  the 
solution  to  large  sysems  of  linear  equations  and  that  of  differential 
equations. 

7.  Rao,  C.  Radhakrishna.  The  Inefficiency  of  Least  Squares:  Extensions  of 
Kantorovich  Inequality.  Technical  Report  #84-47,  Center  for  Multivariate 
Analysis,  November  1984. 

Four  different  measures  of  inefficiency  of  the  simple  least  squares  estimator 
in  the  general  Gauss-Markoff  model  are  considered.  Previous  work  on  the 
bounds  to  some  of  these  measures  is  briefly  reviewed  and  new  bounds  are 
obtained  for  a  particular  measure. 

8.  Bai.  Z.  D..  Krishnaiah,  P.  R.,  and  Yin,  Y.  Q.  On  Limiting  Spectral 
Distribution  of  Product  of  Two  Random  Matrices  when  the  Underlying 
Distribution  is  Isotropic.  Technical  Report  #84-48,  Center  for  Multivariate 
Analysis.  November  1984. 

Let  S  =  XX'  be  distributed  independent  of  a  nonnegative  symmetric  random 
matrix  T,  where  X  =  [xj,  •  •  •  .XjJ:  p  x  n  and  xt.  •  ,xn  is  a  sample  from  an 
isotropic  population  and  the  second  moments  of  the  norm  Ix^l,  (i  = 
1.2.  ■  ■  ■  ,n).  exist.  In  this  paper,  the  authors  proved  that  the  Innit  of  the 
spectral  distribution  of  ST/n  exists. 

9.  Bai.  Z.  D.  A  Note  on  Asymptotic  Joint  Distribution  of  the  Eigenvalues  of  a 
Noncentral  Multivariate  F  Matrix.  Technical  Report  #84-49.  Center  for 
Multivariate  Analysis.  November  1984. 

In  P.  L.  Hsu  (1941 ).  the  proof  of  the  basic  Lemma  3  is  based  on  Lemma  1 
which  is  wrong.  The  aim  of  this  note  is  to  correct  the  proof  of  Lemma  3. 
consequently,  to  ensure  the  main  theorem  in  P.  L.  Hsu  (1941). 

10.  Yin.  Y.  O  Limiting  Spectral  Distribution  for  a  Class  of  Random  Matrices. 
Technical  Report  =<s4-50.  Center  for  Multivariate  Analysis.  December 
1984. 


Let  X  =  [X, 


1.  2.  ...)  be  an  infinite  dimensional  random  matrix.  Tr 


be  a  p  x  p  nonnegative  definite  random  matrix  independent  of  X,  for  p  =  1. 
2,  ...  .  Suppose  -  tr  Tpk— »  H^a.s.  as  p— for  k  =  1,  2,  ...,  and  XH21c‘1/2k 
=oo  .  Then  the  spectral  distribution  of 


where  Xp  =  [X-  :  i  =  1 . p;  j  =  l....n],  tends  to  a  nonrandom  limit 

distribution  as  p  — >  ,°°.n  — » but  p/n  — »  y  >  0,  under  the  mild  conditions 
that  Xy  are  i.i.d.  and  EXj]2  > 

11.  Stoffer.  D.  S.  BoostrapDing  the  Kalman  Filter.  Technical  Report  #84-51, 
Center  for  Multivariate  Analysis,  December  1984. 

The  bootstrap  is  proposed  as  a  method  for  estimating  the  precision  of 
forecasts  and  estimates  of  parameters  of  the  Kalman  Filter  model.  It  is 
shown  that  when  the  system  and  the  filter  is  in  steady  state  the  bootstrap 
applied  to  the  Gaussian  innovations  yields  asymptotically  consistent 
standard  errors.  That  the  bootstrap  works  well  with  moderate  sample  sizes 
and  supplies  robustness  against  departures  from  normality  is  substantiated 
by  empirical  evidence. 

12.  Burbea.  Jacob.  Informative  Geometry  of  Probability  Spaces.  Technical 
Report  #84-  52.  Center  for  Multivariate  Analysis.  December  1984. 

The  paper  is  concerned  with  the  geometrical  properties  that  are  induced  by 
the  local  information  contents  and  structures  of  the  parameter  space  of 
probability  distributions.  Of  particular  interest  in  this  investigation  is  the 
Rao  distance  which  is  the  geodesic  distance  induced  by  the  differential 
metric  associated  with  the  Fisher  information  matrix  of  the  parameter 
space.  Moreover,  following  Efron.  Dawid  and  Amari,  some  affine 
connections  are  introduced  into  the  informative  geometry  of  parameter 
space  and  thereby  elucidating  the  role  of  the  curvature  in  statistical  studies. 
In  addition,  closed  form  expressions  of  the  Rao  distances  for  certain 
families  of  probability  distributions  are  given  and  discussed. 

13.  Bai,  Z.  D..  Knshnaiah,  P.  R..  and  Liang.  W.  Q.  On  Asymptotic  Joint 
Distribution  of  the  Eigenvalues  of  the  Noncentral  MANOVA  Matrix  for 
Nonnormal  Populations.  Technical  Report  #84-53,  Center  for  Multivariate 
Analysis.  December  1984. 

In  this  paper,  the  authors  derived  asymptotic  joint  distribution  of  the 
eigenvalues  of  the  multivariate  analysis  of  variance  matrix  in  the  noncentral 
case  when  the  underlying  distribution  is  not  multivariate  normal. 

14.  Giri.  .V.  and  Sinha.  B.  K.  Robust  Tests  of  Mean  Vector  in  Symmetrical 
Multivariate  Distributions.  Technical  Report  #85-01.  Center  for 
Multiva  :.:'c  Nnaivsis.  January  1985. 

LetY=(Aj  .  •  V  '.  ..Y„')'„VJ'-(.V,/ . A'  )  be  a  n  x  p  random  matrix  with 

probability  demits  function 


where  x  €  x  =  {x:nxp  matrix  I  rank  of  x  =  p } ,  y^Qij,  •  ■  ■  ,\xp)'  e  RP  g  = 
(1,  •  ,1)'  n- vector  and  X  >  0  (positive  definite).  Assume  that  q  e  Q  =  f  q: 

[0,“=),  — *  (0,°°),convex}  and  n  >  p  so  that  X'X  >  0  with  probability  one.  It  is 
proved  that  for  testing  H0:  y.  =  Q  versus  the  alternative  Hj:  y.  #  Q  ,  the 

Hotelling's  T2-test  is  locally  minimax  and  for  testing  Hq:  y^  =  0  versus  the 
alternative  Hp  j^j.  *  0  .  the  appropriate  Hotelling's  T2-test  is  both  UMPI 
and  locally  minimax.  In  the  second  case  y^1p(ji1,...,p^  )',  pj  <  p,  and 

(y.p  +1 . (ip),  X  are  unknown.  The  above  results  generalize  those  of  Giri 
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and  Kiefer  (Ann.  Math.  Statist.,  1964)  under  the  assumption  q 
As  a  technical  tool.  Wijsman's  representation  theorem  is  used. 


normal. 


15.  Bueno  Vanderlei  C.  A  Model  in  Negative  Dependence  Using  Stochastic 
Order.  Technical  Report  #85-02,  Center  for  Multivariate  Analysis, 
February  1985. 

A  model  gives  rise  to  the  NDS  conditions  from  Block,  Savits  and  Shaked. 
The  main  theorem  is:  Let  Xj,  •  •  •  ,  Xjj  be  independent  and  identically 
distributed  random  variables  with  continuous  distribution  F.  Let 

(Fi,....y„)=f(x1....  jc„)\x(k)zz] 

Then  (Y|,  •  •  •  ,  Yn)  is  Negative  Dependent  through  Stochastic  Order 
(NDS). 

The  NDS  condition  implies  the  NUOD  and  NLOD  condition  of  negative 
dependence  that  enable  us  to  get  bounds  in  joint  probabilities. 

The  conditioning  event  considered  occurs  naturally  in  Reliability  Theory  as 
the  time  to  system  failure  for  k-out-of  n  system. 

The  conditioning  in  {g(X(k))<|l}  where  g  is  increasing  and  left  continuous 
function,  also  define  a  NDS  vector.  However  conditioning  in  a  general 
coherent  system  doesn't  rise  a  NDS  vector. 

16.  Khattree.  Ravindra.  and  Yin.  Y.  Q.  An  Inqualitv  and  Its  Application  to  the 
Truncated  Distributions.  Technical  Report  #85-03,  Center  for  Multivariate 
Analysis.  February  1985. 

An  inequality  is  proved  and  its  interpretation  is  given.  Using  the  inequality, 
it  is  shown,  under  some  mild  conditions,  that  for  the  univariate  truncated 
distributions,  the  variance  of  the  truncated  distribution  increases  with  the 
value  of  the  truncation  point. 

17.  Das.  Rita  and  Sinha.  Bimat  K.  Detection  of  Multivariate  Outliers  with 
Dispersion  Slippage  in  F.llipticallv  Symmetric  Distributions.  Technical 
Report  -S5-04.  Center  for  Multivariat  Analysis.  March  1985. 

An  e\ie:t'  -u  .>;  Ferguson's  (Fourth  Berkeley  Symposium  on  Probability 
and  Matrivti-.ut.cul  Stan  -.tics.  1961.  Volume  1)  univariate  normal  results  for 
detection  <•:  .•  it.ier,  with  variance  slippage  is  made  to  the  multivariate 
ellipticniN  ■> : i  i •  ; -:rt .  case  with  dispersion  slippage.  The  locally  optimum 
test  statist:,  we  d-me  possesses  all  three  robustness  properties:  optimality, 
null  and  nonrv.d!  -  Kama  and  Sinha.  Ann.  Statist..  1985).  As  a  technical 
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tool,  Wijsman's  representation  theorem  (Fifth  Berkeley  Symposium  on 
Probability  and  Mathematical  Statistics,  1967,  Volume  1)  is  used. 

18.  Rao,  C.  Radhakrishna.  Evolution  of  Data  Collection  Censuses,  Sample 
Surveys,  and  Design  of  Experiments.  Technical  Report  #85-05,  Center  for 
Multivariate  Analysis,  March  1985. 

Statistics  as  we  understand  and  practice  it  today  has  a  short  history  but  a 
long  antiquity.  Originally  the  word  statistics  meant  data  or  information 
relating  to  a  State,  and  gradually  it  acquired  a  wider  meaning  to  cover 
numerical  data  in  any  field  of  enquiry  and  also  analyses  performed  on  data 
to  draw  inferences  and  to  take  decisions.  In  early  days,  attempts  were  made 
to  make  a  complete  enumeration  when  information  on  people  and  resources 
of  a  country  was  required.  This  was  called  Census.  But  with  increase  in 
demand  for  extensive  and  timely  data  by  governments,  individuals  and 
private  enterprises,  the  method  of  sample  surveys  involving  the  collection 
of  data  only  on  a  small  subset  of  the  population  and  making  estimates  for 
the  whole  population  was  developed.  This  introduced  flexibility,  speed  and 
economy  in  acquiring  information  of  sufficient  accuracy  for  practical 
purposes.  Simultaneously,  methods  of  generating  data  through  well 
designed  experiments  was  developed  with  the  twin  objectives  of  providing 
valid  estimates  of  differences  between  treatments  with  maximum  precision 
for  given  cost.  Design  of  experiments  also  introduced  great  simplicity  in 
the  analysis  and  interpretation  of  data.  Although  census,  sample  surveys 
and  design  of  experiments  were  introduced  for  separate  purposes,  now-a- 
davs  investigations  in  all  fields  of  enquiry  use  a  combination  of  the 
principles  of  all  these  three  methods  of  data  collection.  New  developments 
in  data  collection  are  computer  assisted  telephone  enquiry  (CATI)  and  use 
of  satellites  in  taking  pictures  of  earth  for  estimating  average  and  yield  of 
crops  through  pattern  recognition  techniques. 

19.  Mathew,  Thomas.  On  the  Characterization  of  Nonnegatively  Estimable 
Linear  Combinations  of  Variance  Components.  Technical  Report  #85-06. 
Center  for  Multivariate  Analysis,  March  1985. 

It  is  shown  that,  by  a  reparametrization,  the  problem  of  estimating  a  linear 
combination  of  variance  components  can  be  reduced  to  that  of  estimating  a 
single  variance  component.  Such  a  reduction  is  used  to  obtain  some 
characterizations  of  nonnegatively  estimable  linear  combinations  of 
variance  components.  Characterization  of  nonnegative  estimability  using 
M1NQUE  is  also  discussed. 

20.  Langberg.  Naftali  A.,  Stoffer.  David  S.  Moving  Average  Models  with 
Bivariate  Exponential  and  Geometric  Distributions.  Technical  Report 
#85-07,  Center  for  Multivariate  Analysis.  March  1985. 

Two  classes  of  finite  and  infinite  moving  average  sequences  of  bivariate 
random  vectors  are  considered.  The  first  class  has  bivariate  exponential 
marginals  while  the  second  class  has  bivariate  geometric  marginals.  The 
theory  ot  positive  dependence  is  used  to  show  that  in  various  cases  the  two 
classes  consist  of  associated  random  variables.  Association  is  then  applied 
to  establish  moment  inequalities  and  to  obtain  approximations  to  some  joint 
probabilities  of  the  related  bivariate  point  processes. 
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21.  Rao,  C.  Radhakrishna.  Differential  Metrics  in  Probability  Spaces  Based  on 
Entropy  and  Divergence  Measures.  Technical  Report  #85-08,  Center  for 
Multivariate  Analysis,  April  1985. 

In  this  paper  we  discuss  some  general  methods  of  metrizing  probability 
spaces  through  the  introduction  of  a  quadratic  differential  metric  in  the 
parameter  manifold  of  a  set  of  probability  distributions.  These  methods 
extend  the  investigation  made  in  Rao  (1945)  where  the  Fisher  information 
matrix  was  used  to  construct  the  metric,  and  the  geodesic  distance  was 
suggested  as  a  mesure  of  dissimilarity  between  probability  distributions. 

The  basic  approach  in  the  present  paper  is  first  to  construct  a  divergence  or 
a  dissimilarity  measure  between  any  two  probability  distributions,  and  use 
it  to  derive  a  differential  metric  by  considering  two  distributions  whose 
characterizing  parameters  are  close  to  each  other.  One  measure  of 
divergence  considered  is  the  Jensen  difference  based  on  an  entropy 
functional  as  defined  in  Rao  (1982a).  Another  is  the  f-  divergence  meausre 
studied  by  Csiszar  (1967).  The  latter  class  leads  to  the  differential  metric 
based  on  the  Fisher  information  matrix.  The  geodesic  distances  based  on 
this  metric  computed  by  various  authors  are  listed. 

22.  Krishnaiah.  P.  R..  and  Sarkar.  Shakuntala.  Principal  Component  Analysis 
Under  Correlated  Multivariate  Regression  Equations  Model.  Technical 
Report  #85-09,  Center  for  Multivariate  Analysis,  April  1985. 

In  this  paper,  the  authors  consider  the  problem  of  testing  for  the  equality  of 
the  last  few  eigenvalues  of  the  covariance  matrix  under  correlated 
multivariate  regression  equations  (CMRE)  model.  Asymptotic  distributions 
of  various  test  statistics  are  derived  when  the  underlying  distribution  is 
multivariate  normal.  Some  of  the  distribution  theory  is  extended  to  the  case 
when  the  underlying  distribution  is  elliptically  symmetric. 

23.  Lau.  Ka-Sing,  and  Gu,  Hua-Min.  A  Note  on  an  Integrated  Cauchy 
Functional  Equation.  Technical  Report  #85-10,  Center  for  Multivariate 
Analysis,  April  1985. 

In  characterizing  the  semistable  law,  Shimizu  reduced  the  problem  into 
solving  the  equation 

H(x)=j  H(.x+y >ri(p.-v)(y),  .r>0 

where  p  and  v  are  given  positive  measures  on  10.°°).  In  this  note,  we  obtain 
a  simple  proof  and  show  that  some  of  his  conditions  can  be  weakened. 

24.  Alzaid.  Abdulhamid  A..  Rao.  C.  Radhakrishna  and  Shanbhag.  D.  N.  An 
Extension  of  Spit/er's  Integral  Representation  Theorem  with  ;m 
Application.  Technical  Report  #85-1 1.  Center  for  Multivariate  Analysis. 

Using  a  ivj'.s  approach,  an  extended  version  of  Spit/.er's  integral 
representation  for  stationary  measures  of  a  discrete  branching  process  is 
obtained  1:  ; ,  remit  in  used  to  provide  a  complete  solution  to  a  problem  in 
damage  modeA  mtisfung  a  generalized  Rao-Rubin  condition. 

25.  Singh.  Radhe\  S  and  Ahmad.  Manzoor.  Modified  Nonparametric  Kernel 


Estimates  of  a  Regression  Function  and  Then  Consistencies  with  Rates. 
Technical  Report  #85-12.  Center  for  Mulnvanate  Analysis.  April  1985. 

Two  sets  of  modified  kernel  estimates  of  a  regression  function  are 
proposed:  one  when  a  bound  on  the  regression  funcnon  is  known  and  the 
other  when  nothing  of  this  sort  is  at  hand.  Explicit  bounds  on  the  mean 
square  errors  of  the  estimators  are  obtained.  Pointwise  as  well  as  uniform 
consistency  in  mean  square  and  consistency  in  probability  of  the  estimators 
are  proved.  Speed  of  convergence  in  each  case  is  invesngated. 

26.  Alzaid,  Abdulhamid  A..  Rao,  C.  Radhaknshna  and  Shanbhag,  D.  N.  An 
Application  of  the  Perron-Frobenius  Theorem  to  a  Damage  Model  Problem. 
Technical  Report  #85-13.  Center  for  Mulnvanate  Analysis.  Apnl  1985. 

Using  the  Perron-Frobenius  theorem,  it  is  established  that  if  (X.Y)  is  a 
random  vector  of  non-negative  integer  valued  components  such  that  Y  <  X 
almost  surely  and  two  modified  Rao-Rubm  conditions  hold,  then  under 
some  mild  assumptions  the  distribution  of  (X.Y)  is  uniquely  determined  by 
the  conditional  distribution  of  Y  given  X.  This  result  extends  the  recent 
unpublished  work  of  Shanbhag  and  Taillie  ( 1979)  on  damage  models. 

27.  Krishnaiah.  P.  R..  Lin.  J.  and  Wang.  L.  Inference  on  the  Ranks  of  the 
Canonical  Correlation  Matrices  for  Elliptically  Symmetric  Populations. 
Technical  Report  #85-  14,  Center  for  Multivariate  Analysis,  April  1985. 

In  this  paper,  the  authors  considered  the  likelihood  ratio  tests  and  some 
other  tests  for  the  ranks  of  the  canonical  correlation  matrices  when  the 
underlying  distributions  are  reaal  and  complex  elliptically  symmetric 
distributions.  Also  asymptotic  joint  distributions  of  the  eigenvalues  of  the 
sample  canonical  correlation  matrices  are  derived  under  the  assumptions 
mentioned  above  regarding  the  underlying  distributions.  Finally, 
applications  of  tests  for  the  rank  of  the  complex  canonical  correlation 
matrix  in  the  area  of  time  series  in  the  frequency  domain  are  discussed. 

28.  Bai,  Z.  D.  Exponential  Bound  for  Error  Probability  in  NN-Discrimination. 
Technical  Report  #85-15,  Center  for  Multivariate  Analysis.  April  1985. 

Let  (0j.Xj).  •  •  •  ,(0n,Xn)  be  n  simple  samples  drawn  from  the  population 

(0.X)  which  is  a  (  1.2 . s}  x  Revalued  random  vector.  Suppose  that  X  is 

known.  Let  Xn'  be  the  nearest  one  among  A', _ Xr,  from  X,  in  the  sense  of 

Euclidean  norm  or  norm,  and  let  0n  be  the  0-value  paired  with  Xn\  The 
posterior  error  probability  is  defined  by  Ln=P(Qn'*  01 
{(0j.Xj),  ■  .( 0n.Xn) } .  It  is  well  known  that  E L,=P(Qn'  *  0)  has  always  a 

limit  R.  In  this  paper  it  is  shown  that  for  any  e  >  0  there  exists  constant  C  < 
°°  and  b  >  0  such  that  P(  I  Lr-R  l>£)  <  ce'^n.  under  the  only  assumption 
that  tire  marginal  distribution  of  X  is  nonatomic. 

29.  Al/.aid.  Abdulhamid  A..  Rao.  C.  Radhaknshna.  and  Shanbhag. 
D  N.  Characterization  of  Discrete  Probability  Distributions  by  Partial 
Independence.  Technical  Report  #85-16.  Center  for  Multivariate  Analvsis, 
May  1985. 

If  X  and  Y  are  random  variables  such  that  P(X  >  Y')  =  1  and  the  conditional 
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distribution  of  Y  given  X  is  binomial,  then  Moran  (1952)  showed  that  Y 
and  (X-Y)  are  independent  if  X  is  Poisson.  We  extend  Moran’s  result  to  a 
more  general  type  of  conditional  distribution  of  Y  given  X,  using  only 
partial  independence  of  Y  and  X-Y.  This  provides  a  generalization  of  a 
recent  result  of  Janardhan  and  Rao  (1982)  on  the  characterization  of 
generalized  Polya- Eggenberger  distribution.  A  variant  of  Moran's  theorem 
is  proved  which  generalizes  the  results  of  Patil  and  Seshadri  ( 1 964)  on  the 
characterization  of  the  distribution  of  a  random  variable  X  based  on  some 
conditions  on  the  conditional  distribution  of  Y  given  X  and  the 
independence  of  Y  and  X-Y. 

30.  Yin,  Y.  Q.  and  Bai.  Z.  D.  Spectra  for  Large  Dimensional  Random  Matrices. 
Technical  Report  #85-17,  Center  for  Multivariate  Analysis,  May  1985. 

In  this  paper,  the  authors  reviewed  some  recent  developments  in  the  area  of 
large  dimensional  random  matrices. 

31.  Hanumara,  R.  Choudary.  An  Alternate  Derivation  of  the  Distribution  of  the 
Conditioned  Signal-To-Noise  Ratio.  Technical  Report  #85-18,  Center  for 
Multivariate  Analysis,  May  1985. 

A  result  in  multiple  regression  analysis  is  used  to  derive  the  probability 
distribution  of  the  conditioned  signal-to-noise  ratio. 

32.  Alzaid.  Abdulhamid  A.,  Rao,  C.  Radhakrishna  and  Shanbhag,  D.  N.  A  Note 
on  the  Integrated  Cauchy  Functional  Equation.  Technical  Report  #85-19, 
Center  for  Multivariate  Analysis.  May  1985. 

New  proof  is  given  using  de  Finneti's  theorem  on  infinite  sequence  of 
exchangeable  random  ariables  for  the  solution  of  the  integrated  Cauchy 
functional  equation  studied  by  Lau  and  Rao  (1982)  in  the  continuous  case 
and  by  Shanbhag  (1977)  in  the  discrete  case. 

33.  Fang,  C.,  and  Krishnaiah,  P.  R.  On  Asymptotic  Distribution  of  the  Test 
Statistic  for  the  Mean  of  the  Non-Isotropic  Principal  Component. 
Technical  Report  #85-20,  Center  for  Multivariate  Analysis,,  May  1985. 

In  this  paper,  the  authors  derived  the  large  sample  distribution  of  the  t 
statistic  based  upon  the  observations  on  the  first  principal  component 
instead  of  the  original  variables.  It  is  shown  that  he  above  statistic  is 
distributed  asymptotically  as  Student's  t  distribution. 

34.  Ahmad,  Manzoor,  and  Singh,  Radhev  S.  Role  of  Axillary  Variate  and 
Additional  Data  in  Density  Estimation.  Technical  Report  #85-21,  Center 
for  Multivariate  Analysis,  May  1985. 

Some  new  estimators  of  a  univariate  p.d.f.  by)  of  a  random  variable  Y. 
based  on  a  set  of  observations  taken  from  a  bivariate  joint  densitv  6(x.y )  of 
Y  and  a  suitably  chosen  concomitant  variable  X.  have  been  investigated. 
Asymptotic  unbiasedness,  mean  square  consistency,  asymptotic  normality 
and  rates  of  convergence  have  been  established.  A  related  problem  of 
estimation  of  a  conditional  density  has  also  been  studied. 

35.  Alzaid.  Abdulhamid  A..  Rao.  C.  Radhaknshna.  and  Shanbhag.  D.  N.  A 
Variant  of  Shanbhag 's  Lemma  Arising  out  of  a  Modified  Rao-Rubin 
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Condition.  Technical  Report  #85-22.  Center  for  Multivariate  Analysis, 
May  1985. 

A  modified  Rao-Rubin  condition  for  damage  models  gives  rise  to  a 
recurrence  relation  which  is  somewhat  different  from  that  considered  by 
Shanbhag  (1977).  A  complete  solution  to  the  new  recurrence  relation  is 
obtained  and  its  applications  are  indicated. 

36.  Ahmad  Manzoor.  New  Measures  of  Diversity.  Technical  Report  #85-23. 
Center  for  Multivariate  Analysis.  June  1985. 

The  problem  of  measuring  diversity  within  populations  and  dissimilarity  or 
similarity  between  populanons  has  been  extensively  treated  in  the  literature. 
In  this  context  a  general  procedure  called  Analysis  of  Diversity  has  been 
outlined  and  examined  by  C.  R.  Rao  in  a  series  of  papers. 

In  this  paper  we  propose  three  new  measures  of  diversity  and  study  related 
inference  problems.  Denote  by  Sk  the  simplex  Sk  =(&&s(7iI,  ■  ■  ■  ,nkY, 
7tj>0.  Dtj  =  1  }.  Then  the  proposed  measures  are  of  the  form: 

1  -amI;Tt;0,„( 7ty).  m=1.2.3  where  Oj(x)  =  (l+k'1-x)‘7,  y  >0,  p,(x)  = 
(2-k'7-x‘0*1.y>  0,  0j(x)  =  (a?+(  l-x)7)'',  0  <  y<  1.  and  the  a’s  are  suitable 
normalizing  constants.  Estimation  of  H^jp.  derivation  of  the  penalty 
function  and  cross  entropy  and  the  problem  of  testing  independence  have 
been  treated.  Asymptotic  distributions  of  relevant  test  statistics  are 
indicated. 

37.  Ahmad,  Manzoor.  and  Singh.  R.  S.  Double  Sampling  in  Estimation  of  A 
Ratio.  Technical  Report  #85-24.  Center  for  Multivariate  Analysis.  June 
1985. 

Problems  connected  with  the  estimation  of  the  ratio  of  the  means  of  a  finite 
bivariate  population  have  been  considered  in  this  report.  The  usual 
estimator  of  the  ratio,  based  on  the  means  of  a  bivariate  simple  random 
sample  drawn  without  replacement  (s.r.s.  (w.o.r.)).  has  been  compared  with 
estimators  based  on  alternative  double  sampling  design.  Under  this  design 
a  very  large  s.r.s.  (w.o.r.)  is  drawn  for  measunng  only  one  of  the  variables 
and  a  subsample  (s.r.r.f w.o.r.))  is  drawn  out  of  the  first  phase  units  for 
measunng  the  other  variable.  Efficiency  and  bias  compansons  have  been 
made  by  subjecting  each  of  the  competitors  to  the  same  budgetary 
constraint.  It  turns  out  that  deviation  from  the  usual  set-up  sometimes  leads 
to  better  sampling  'strategies'. 

38.  Mathew.  Thomas.  Hermitian  and  Nonnegativity  Preserving  Subspaces. 
Technical  Report  =*85-25.  Center  for  Multivariate  Analysis.  July  1985. 

Tie  concept'  of  hermitian  preserving  and  nonnegalivity  preserving 
subspaces  of  complex  square  matrices  are  introduced.  Characterizations  of 
such  sub-paces  are  obtained  and  applications  are  discussed. 

39.  Baiqi.  Mian,  and  Limheng.  Zhao.  Non-uniform  Bounds  of  Normal 
Approximation  lor  bmite-population  U-statisncs.  Technical  Report  #85-26. 
Center  for  Mult. variate  Analysis.  July  1985. 
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Let  An  be  a  population  with  N  balls  bearing  numbers  aN1,  •  •  • , 
respectively.  Draw  n  balls  from  AN  randomly  without  replacement,  and 
denote  the  numbers  appearing  on  these  n  balls  by  Xj,  ••  •  ,Xn.  Suppose  that 
<j)N(x,y)  be  a  Borel-measurable  function,  symmetric  in  x  and  y.  Set 

^  =  I  ^XjJ(k),QN=E^X  ,Jf2). 

1  £)<k<,n 

g(X1)  =  £t<D^1^2)IXI) 

ag“=Var{g(X[)l.  In  this  paper  we  established  that,  if  there  exists  fixed 
constants  Xl  and  such  that  0  <  >,,  <  n/N  <  h,  <  1.  then  it  is  valid  for  all 
positive  integer  n  and  real  x  that 

mun-eN) 

!/>(— <xH>(x)l  <  Cn~m(5  -3£I<MX,.X,)I3(1+Lxir3 
2\1  -n/Nog 

where  <&(x)  is  the  standard  normal  distribution  function,  and  C  is  an 
absolute  constant  depending  solely  on  ^  and 

40.  Bimal  K.  Sinha  and  Malay  Ghosh.  Inadmissibility  of  the  Best  Equivariant 
Estimators  of  the  Variance-Covanance  Matrix  and  the  Generalized 
Vanance  Under  Entropy  Loss  Technical  Report  No.  85-27,  Center  for 
Multivariate  Analysis,  July  1985. 

Based  on  a  data  matrix  X  =  (X  j.  •  ,X  k):  px  k  with  independent 
columns 

X  |  ~  Np(4  j.  .£)•  and  an  mdependent  Wishart  matrix  S:  pxp  ~  W  (n.Z), 
estimators  dominating  the  beet  equivariant  estimators  of  I  and  El  are 
obtained  under  two  types  of  entropy  loss.  For  simultaneous  estimation  of 
the  mean  vector  and  the  variance  covariance  matrix  of  a  multinormal 
population,  a  suitable  entropy  loss  is  developed  and  testimators  dominating 
the  pair  consisting  of  the  sample  mean  vector  and  the  best  multiple  of  the 
sample  Wishart  matrix  are  derived.  A  technique  of  Sinha  (Jour.  Mult. 
Analysis,  1976)  is  heavily  exploited. 

41.  Demetrios  G.  Kaffes,  M.  Bhaskara  Rao.  Thomas  Mathew  and 
K.  Subramanyam.  On  The  Matrix  Convexity  of  the  Moore- Penrose  Inverse 
and  Some  Applications.  Technical  Report  No.  85-28.  Center  for 
Multivariate  Analysis.  July  1985. 

It  is  well  known  that  if  A  and  B  are  two  positive  definite  matrices  of  the 
same  order  and  ()<  >.<  1,  then 

1/dt  i  1— >.u5|  <>„4"‘+(  1  -k)R  1 

It  in  easy  to  construct  an  example  consisting  of  two  positive  semi-definite 
matrices  tor  \\  inch  the  above  inequality 

IaA-m  1  l+</_4++(l-/.)£+ 

for  every  ()</.<  1 .  As  an  application,  we  give  a  sufficient  condition  under 
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which  the  inequality  (EA)+  <  E(A+)  is  valid,  where  A  is  a  square  matrix  of 
random  variables  which  is  almost  surely  positive  semi-definite, 
generalizing  the  the  well-known  result  (EA)'1  <  EA'1  when  A  is  almost 
surely  positive  definit. 

42.  Bai,  Z.  D.,  and  Liang,  W.  Q.  Strong  Representation  of  Weak  Convergence. 
Technical  Report  No.  85-29.  Center  for  Multivariate  Analysis,  July  1985. 

Let  (in.  n  =  1,2 . and  p.  be  a  given  sequence  of  probability  measures  each 

of  which  is  defined  on  a  complete  separable  metric  space  Sn  and  S 
respectively.  Also,  a  sequence  of  measurable  mappings  qn  from  Sn  into  S  is 
given.  In  this  paper,  it  is  proved  that  if  pn  o  (jf1  weakly  converge  to  p,  then 

there  is  a  probability  space  (Q.F.P),  on  which  we  can  define  a  sequence  of 
random  elements  Xn,  from  into  Sn.  and  a  random  element  X,  from  Q  into 
S,  such  that  p^  is  the  distribution  of  X^,  p  is  the  distribution  of  X  and 
lim  db^X^)  =  X  pointwise.  The  result  of  Skorokhod  (1956)  is  a  special 

case  of  the  result  of  this  paper.  Some  applications  in  the  area  of  random 
matrices,  etc.,  are  also  given. 

43.  Lin-Cheng  Zhao.  An  Inequality  Concerning  the  Deviation  Between 
Theoretical  and  Empirical  Distributions.  Technical  Report  No.  85-30, 
Center  for  Multivariate  Ana  Analysis,  August  1985. 

In  this  paper  the  author  established  an  inequality  concerning  the  uniform 
deviation  between  theoretical  and  empirical  distributions.  An  application  in 
strong  convergence  of  nearest  neighbor  density  estimate  is  also  discussed. 

44.  Fujikoshi,  Y.,  Krishnaiah,  P.  R.,  and  Schmidhammer,  J.  Effect  of 
Additional  Variables  in  Principal  Component  Analysis.  Discriminant 
Analysis  and  Canonical  Correlation  Analysis.  Technical  Report  No.  85-31, 
Center  for  Multivariate  Analysis,  August  1985. 

In  this  paper,  the  authors  derived  asymptotic  distributions  of  changes  in 
certain  functions  of  the  eigenvalues  of  the  sample  covariance  matrix. 
MANOVA  matrix  and  canonical  correlation  matrix  when  some  variables 
are  added  to  the  original  sets  of  variables.  The  above  results  are  useful  in 
finding  out  as  to  whether  the  new  variables  give  additional  information. 

45.  Takeaki  Kariya.  ;uid  Bimal  K.  Sinha.  Optimality  Robustness  of  Tests  in 
Two  population  Problems.  Technical  Report  No.  85-32.  Center  For 
Multivariate  Analysis.  August  1985. 

The  problems  of  testing  the  equality  of  two  location  parameters  (without 
the  presence  of  scale  parameters)  and  the  equality  of  two  scale  parameters 
(with  or  orient  the  presence  of  location  parameters)  based  on  two 
independent  '.unples  are  considered.  It  is  shown  that,  under  some  mild 
condition-.  .  me  underlying  distributions,  the  standard  optimum  invariant 
tests  de:  .  1  ..-  ier  normality  or  exponentially  continue  to  be  optimum  for  a 
wide  class  . •:  distributions.  However,  the  tests  are  not  null  robust.  As  a 
technical  tool.  Wiisman's  (Fifth  Berkeley  Symposium.  Vol.  1.  1967) 
representation  theorem  is  used. 


17 


46.  Kariya,  T,  Fujikoshi,  Y.  and  Krishnaiah,  P.  R.  On  Tests  for  Selection  of 
Variables  and  Independence  Under  Multivariate  Regression  Model. 
Technical  Report  No.  85-  33,  Center  for  Multivariate  Analysis,  August 
1985. 

In  this  paper,  the  authors  consider  various  procedures  for  testing  the 
hypotheses  of  independence  of  two  sets  of  variables  and  certain  regression 
coefficients  are  zero  under  the  classical  multivariate  regression  model. 
Various  properties  of  these  procedures  and  the  asymptotic  distributions 
associated  with  these  procedures  are  also  considered. 

47.  Rao,  Radhakrishna  C.,  and  Rao  C.  Veerendra.  Computation  of  the 
Stationary  Values  of  the  Product  of  Two  Raleigh  Quotients.  Technical 
Report  No.  85-34,  Center  for  Multivariate  Analysis.  September  1985. 

A  computational  algorithm  is  developed  for  finding  the  stationary  values  of 
the  function  x’CxKx'  Ax) 1  ^(x'Bx) 1  ^  where  A  and  B  are  positive  definite  and 
C  is  is  a  symmetric  matrix .  The  square  of  the  function  under  consideration 
is  the  product  of  two  Raleigh  coefficients  xCx/x'Ax  and  x'Cx/xBx.  The 
general  problem  occurs  in  multivariate  analysis  in  the  computation  of 
homologous  canonical  variates  in  studying  relationships  between  two  sets 
of  homologous  measurements.  The  special  case  with  C  =  I  occurs  in 
designing  control  systems  with  minimum  norm  feedback  matrices. 

48.  Rosenblatt-Roth,  M.  The  Relative  Entropy  of  a  Random  Vector  with 
Respect  to  Another  Random  Vector.  Technical  Report  No.  85-35,  Center 
for  Mulnvariate  Analysis  October  1985. 

In  this  paper,  the  author  discusses  some  problems  connected  with  the 
concept  of  entropy.  Part  of  the  paper  is  expository  in  nature  and  the 
remaining  material  deals  with  new  results  obtained  by  the  author. 

49.  Krishnaiah.  P.  R  ,  Lin,  J.,  and  Wang,  L.  Tests  for  the  Dimensionality  of  the 
Regression  Matrices  when  the  Underlying  Distributions  are  Elliptically 
Symmetric.  Technical  Report  No.  85-36,  Center  for  Multivariate  Analysis. 
October  1985. 

In  this  paper,  the  authors  derive  likelihood  ratio  tests  for  the  dimensionality 
of  the  regression  matrices  for  the  cases  when  the  jomt  distributions  of  the 
observations  are  real  and  complex  elliptically  symmetric.  The  authors  also 
derive  asymptotic  distributions  of  the  above  test  statistics  for  two  situations. 
In  the  first  situation,  the  joint  distribution  of  the  observations  is  elliptically 
symmetric  whereas  the  second  situation  assumes  that  the  observations  are 
distributed  independently  as  elliptically  symmetric. 

50.  Zhao.  L.  C..  Krishnaiah.  P.  R..  and  Bai.  Z.  D.  On  Detection  of  Number  of 
Signals  in  Presence  of  White  Noise.  Technical  Report  No.  85-37.  Center 
for  Multivariate  Analysis.  October  1985. 

In  this  paper,  the  authors  propose  procedures  for  the  detection  of  the 
number  of  signals  in  presence  of  Gaussian  white  noise.  The  methods  used 
fall  within  the  framework  of  the  model  selection  procedures  using  an 
information  theoretic  criterion  The  strong  consistency  of  the  estimates  of 
the  number  of  signals,  under  different  situations,  is  established.  Extensions 
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of  the  results  are  also  discussed  to  the  case  when  the  noise  is  not  necessarily 
Gaussian. 


51. Khatri,  C.  G.,  and  Rao,  C.  Radhakrishna.  Effects  of  Estimated  Noise 
Covariance  Matrix  in  Optimal  Signal  Detection.  Technical  Report  No. 
85-38.  Center  for  Multivariate  Analysis.  October  1985. 

There  is  loss  of  efficiency  when  an  estimated  noise  covariance  matrix  is 
used  in  the  place  of  the  unknown  true  noise  covariance  matrix  in  the 
construction  of  the  optimum  filter  for  signal  detection.  In  the  case  of 
detecting  a  single  signal  specified  by  a  real  or  a  complex  vector,  we 
investigate  the  extent  of  this  loss  by  obtaining  an  exact  confidence  bound 
for  the  realized  signal  to  noise  ratio.  We  also  give  an  estimate  of  this  rano 
which  is  useful  in  optimum  selection  of  features.  Some  of  these  results  are 
extended  to  the  case  of  discrimination  between  a  number  of  given  signals. 

52.  Zhao.  L.  C.,  and  Krishnaiah.  P.  R.,  and  Bai.  Z.  D.  On  Detection  of  the 
Number  of  Signals  When  the  Noise  Covariance  Matrix  is  Arbitrary. 
Technical  Report  No.  85-  39,  Center  for  Multivariate  Analysis,  October 
1985. 

In  this  paper,  the  authors  proposed  model  selection  methods  for 
determination  of  the  number  of  signals  in  presence  of  noise  with  arbitrary 
covariance  matrix.  This  problem  is  related  to  finding  the  multiplicity  of  the 
smallest  eigenvalue  of  ZsZj1.  where  Z>  =  T  +  A.Z],  E,  and  Z.  are 

covariance  matrices,  /„  is  a  scalar,  and  T  is  non-negative  definite  matrix  and 
is  not  of  full  rank.  Also,  the  authors  proposed  methods  for  determination  of 
the  multiplicities  of  various  eigenvalues  of  Z^Z'1.  The  methods  used  in 

these  procedures  are  based  upon  certain  information  theoretic  criteria.  The 
strong  consistency  of  these  criteria  is  established  in  this  paper. 

53.  Khatri.C.  G..  Rao,  C.  Radhaknshna,  and  Sun.  Y.  N.  Tables  for  Obtaining 
Confidence  Bounds  for  Realized  Signal  to  Noise  Ratio  with  An  Estimated 
Discriminant  Function.  Technical  Report  No.  85-40.  Center  for 
Multivariate  Analysis.  November  1985. 

Percentage  points  of  a  new  distribution  involving  a  confluent- 
hypergeometric  distribution  obtained  by  Khatri  and  Rao  (1985)  are 
tabulated.  The  use  of  the  tabulated  values  in  obtaining  a  lower  confidence 
bound  for  the  realized  signal  to  noise  ratio  based  on  an  estimated 
discriminant  function  for  signal  detection  is  explained. 

54.  Duilami.  N\.  Subramanyam.  K.  and  Rao.  M.  Bhaskara.  Selecting  the  Best 

Population:  A  Decision  Theoretic  Approach:  The  Case  of  Pareto 

Distribution.  Technical  Report  No.  85-41.  Center  for  Multivariate 
Analysis.  November  ll)S5. 

A  decision  theoretic  approach  is  followed  m  selecting  the  best  of  k  Pareto 
population-,  taking  into  account  the  cost  of  sampling  and  penalities  for 
wrong  decisions  Mumnax  sample  sizes  are  determined  under  various 
types  of  penality  functions. 

55.  Dailanu.  N  .  Rao  M.  Bhaskara.  and  Subramanyam.  K.  On  the  Selection  of 


Best  Gamma  Population:  Determination  of  Minimax  Sample  Sizes. 
Technical  Report  No.  85-42,  Center  for  Multivariate  Analysis,  November 
1985. 

Selecting  the  best  Gamma  population  lrom  a  given  set  of  Gamma 
populations  is  treated  from  a  decision  theoretic  point  of  view.  Cost  of 
sampling  and  penalities  for  wrong  decision  play  a  role  in  the  determination 
of  optimum  common  sample  sizes.  Minimax  sample  sizes  are  determined 
under  two  different  penalty  functions. 

56.  Zhao,  L.  C..  Krishnaiah.  P.  R..  and  Bai.  Z.  D.  Remarks  on  Certain  Criteria 
for  Detection  of  Number  of  Signals.  Technical  Report  No.  85-43.  Center 
for  Multivariate  Analysis,  November  1985. 

In  this  note,  we  derive  the  asymptotic  distribution  of  logarithm  of  the 
likelihood  ratio  statistic  for  testing  the  hypothesis  that  the  number  of  signals 
is  equal  to  q  against  the  alternative  that  it  is  equal  to  k  (specified)  for  a 
special  case.  This  distribution  is  not  chi-square.  The  above  statistic  also 
arises  (see  Wax  and  Kailath  (1985))  in  studying  consistency  property  of 
MDL  and  AIC  criteria  for  detection  of  the  number  of  signals. 

57.  Zhao,  L.  C.  Exponential  Bounds  of  Mean  Error  for  the  Nearest  Neighbor 
Estimates  of  Regression  Functions.  Technical  Report  No.  85-44,  Center  for 
Multivariate  Analysis.  November  1985. 

Let  (X.Y).  (Xj.Yj),  ,(Xn,Yn)  be  i.i.d.  Rr  x  R-  valued  random  vectors 
with  EIYI<°°.  and  let  mn(x)  be  a  nearest  neighbor  estimate  of  the  regression 
function  mtx)  =  E(YlX=x).  In  this  paper,  we  establish  an  exponential 
bound  of  the  mean  deviation  between  mn(x)  and  m(x)  given  the  training 

sample  Zn  =  (X,.Yj.  •  •  •  .X^Y^  under  the  conditions  as  weak  as  possible. 
This  is  a  substantial  improvement  on  Beck’s  Result. 

58.  Bai.  Z.  D..  Knshnaiah.  P.  R..  and  Zhao.  L.  C.  On  Rates  of  Convergence  of 
Efficient  Detection  Criteria  in  Signal  Processing  with  With  Noise. 
Technical  Report  No. 85-  45.  Center  for  Multivariate  Analvsis,  November 
1985. 

Zhao.  Krishnaiah  and  Bai[6]  proposed  certain  information  theoretic  criteria 
for  detection  of  the  number  of  signals  under  an  additive  model  with  white 
noise  when  the  noise  variance  is  known  or  unknown.  They  showed  that  the 
above  criteria  are  strongly  consistent  even  when  the  underlying  distribution 
is  not  necessarily  Gaussian.  In  this  paper  the  authors  obtained  upper 
bounds  on  the  probabilities  of  correct  detection. 

59.  Stoffer.  David  S  Walsh-Fourier  Analysis  of  Discrete- Valued  Time  Senes. 
Techrh.i!  Report  \o  X5-46.  Center  for  Multivariate  Analysis.  November 
Lhs5. 

An  appr.'.u::  •  ::.c  .Utah  ses  of  discrete-valued  time  senes  is  discussed. 

The  .me  A  ....  mphshed  in  the  spectral  domain  using  the  Walsh- 

Fourier  transt  irm  which  is  based  on  Walsh  functions.  This  approach  will 
enable  an  m\ estimator  of  discrete  systems  to  analyze  the  data  in  terms  of 
square-w.r.  etorm>  .md  sequence  rather  than  sine-waves  and  frequency.  We 


develop  a  general  signal-plus-noise  type  model  for  discrete-valued  time 
series  in  which  Walsh-Fourier  spectral  analysis  is  of  interest.  We  consider 
the  problems  of  detecting  whether  or  not  a  common  signal  exists  in 
repeated  measures  on  discrete-valued  time  series  and  in  discrete- valued 
processes  collected  in  an  experimental  design.  We  show  that  these  models 
may  depend  on  unknown  regression  parameters  and  we  develop  consistent 
estimates  of  these  parameters  based  on  the  finite  Walsh-Fourier  transform. 
Applications  to  certain  Markov  models  are  given,  however,  the  methods 
presented  also  apply  to  non-Markov  cases. 

60.  Khatri,  C.  G.,  and  Rao,  C.  Radhakrishna.  Test  for  a  Specified  Signal  When 
The  Noise  Covariance  Matrix  is  Unknown.  Technical  Report  No. 85-47, 
Center  for  Multivariate  Analysis,  November  1985. 

In  the  univariate  case  it  is  well  known  that  the  one  sided  t  test  is  uniformly 
most  powerful  for  the  null  hypothesis  against  all  one  sided  alternatives. 
Such  a  property  does  not  easily  extend  to  the  multivariate  case.  In  this 
paper,  a  test  is  derived  for  the  hypothesis  that  the  mean  of  a  vector  random 
variable  is  zero  against  specified  alternatives,  when  the  covariance  matrix  is 
unknown.  The  test  depends  on  the  given  alternatives  and  is  more  powerful 
than  Hotelling's  T2.  The  results  are  derived  both  for  real  and  complex 
vector  observations  and  under  normal  and  spherical  distributions.  The 
properties  of  the  proposed  tests  are  investigated  in  detail  when  a  single 
alternative  is  specified. 

61.  Rao.  M.  Bhaskara,  Krishnaiah,  P.  R„  and  Subramanyam,  K.  A  Structure 
Theorem  on  Bivariate  Positive  Quadrant  Dependent  Distributions  and  Tests 
for  Independence  in  Two-Way  Contingency  Tables.  Technical  Report  No. 
85-48,  Center  for  Multivariate  Analysis,  December  1985. 

In  this  paper,  the  set  of  all  bivariate  positive  quadrant  dependent 
distributions  with  fixed  marginals  is  shown  to  be  compact  and  convex. 
Extreme  points  of  this  convex  set  are  enumerated  in  some  specific 
examples.  Applications  are  given  in  testing  the  hypothesis  of  independence 
against  strict  positive  quadrant  dependence  in  the  context  of  ordinal 
contingency  tables.  Various  procedures  based  upon  certain  functions  of  the 
eigenvalues  of  a  random  matrix  are  also  proposed  for  testing  for 
independence  in  two-way  contingency  table.  The  performance  of  some 
tests  one  of  w  hich  is  based  on  eigenvalues  of  a  random  matrix  is  compared. 

62.  Taiuguchi.  Masanobu.  Validity  of  Edgeworth  Expansions  of  Minimum 
Contrast  Estimators  for  Gaussian  ARMA  Processes.  Technical  Report  No. 
85-49.  Center  for  Multivariate  Analysis.  December  1985. 

Let  { X, }  be  a  Gaussian  ARMA  process  with  spectral  density  fQ( /.)  .  where 
G  is  an  unknown  parameter.  To  estimate  0  we  propose  a  minimum  contrast 
estimation  method  which  includes  the  maximum  likelihood  method  ;uid  the 
quasi-maximum  likelihood  method  as  special  cases.  Let  be  the  minimum 
contrast  estimator  of  0.  ^ken  we  derive  the  Edgeworth  expansion  of  the 
distribution  of  0  t  up  to  third  order,  and  prove  its  validity.  By  this 
Edgeworth  expansion  we  can  see  that  this  minimum  contrast  estimator  is 
always  second-order  asymptotically  efficient  in  the  class  of  second-order 


asymptotically  median  unbiased  estimators.  Also  the  third-order 
asymptotic  comparisons  among  minimum  contrast  estimators  will  be 
discussed. 

63.  Zhao,  L.  C.  Exponential  Bounds  of  Mean  Error  for  the  Kemal  Estimates  of 
Regression  Functions.  Technical  Report  No.  85-50.  Center  for  Multivariate 
Analysis,  December  1985. 

Let  (X,Y).  (Xj,  Y,).  •  •  •  .(Xjj.Yj,)  be  i.i.d.  RrxR-valued  random  vectors  with 
EIYI  <<».  and  let  Qn(x)  be  a  kernel  estimate  of  the  regression  function  Q(x) 
=  E(YIX  =  x).  In  this  paper,  we  establish  an  exponential  bound  of  the  mean 
deviation  between  Qn(x)  and  Q(x)  given  the  training  sample  ZT  =  (Xj.Yj, 

■  •  ■  .Xjj.Yjj),  under  the  conditions  as  weak  as  possible. 

64.  Rao,  M.  Bhaskara,  Krishnaiah.  P.  R..  and  Subramanyam.  K.  Extreme  Point 
Method  in  The  Determination  of  the  Structure  of  a  Class  of  Bivariate 
Distributions  and  Some  Applications  to  Contingency  Tables.  Technical 
Report  No.  86-01,  Center  for  Multivariate  Analysis,  January  1986. 

A  decomposition  of  the  class  of  all  bivariate  positive  quadrant  dependent 
distributions  into  compact  convex  subsets  is  obtained.  Extreme  points  of 
these  compact  convex  subsets  are  investigated.  A  similar  decomposition 
works  out  for  the  class  of  all  bivariate  distributions.  Extreme  points  of  the 
compact  convex  sets  figuring  in  this  decomposition  are  analyzed.  An 
application  to  contingency  tables  is  presented. 

65.  Stoffer,  David  S.  Bootstrapping  Nonlinear  Least  Squares  Estimates  in  the 
Kalman  Filter  Model.  Technical  Report  No.  86-02,  Center  for  Multivariate 
Analysis.  January  1986 

The  bootstrap  is  proposed  as  a  method  for  estimating  the  precision  of 
forecasts  and  maximum  likelihood  estimates  of  the  transition  parameters  of 
the  Kalman  filter  model  when  the  estimates  are  obtained  via  Newton- 
Raphson.  It  is  shown  that  when  the  system  and  the  filter  are  in  steady-state, 
the  bootstrap  applied  to  the  Gaussian  innovations  yields  asymptotically 
consistent  standard  errors.  That  the  bootstrap  works  well  with  moderate 
sample  sizes  and  supplies  robustness  against  departures  from  normality  is 
substantiated  by  empirical  evidence. 

66.  Bai.  Z.  D..  Krishnaiah.  P.  R..  and  Zhao.  L.  C.  Signal  Processing  Using 
Model  Selection  Methods.  Technical  Report  No.  86-03.  Center  for 
Multivariate  Analysis.  January  1986. 

In  this  paper,  the  authors  gave  a  review  of  some  recent  developments  on 
multivariate  statistical  techniques  for  detection  of  the  number  of  signals 
using  eigene  .dues.  The  main  emphasis  of  the  review  is  on  EFFICIENT 
DETE'IK'N  CRITERION  (EDO  procedures  proposed  by  the  authors 
recent !\ .  Uv  e  procedures  are  strongly  consistent. 

67.  Rao.  (.'  R..  "  An -linn  A  NODI  V:  Generalization  of  A  NOVA  Through 

En trope  .mu  Cross  Entropy  Functions.  Technical  Report  No.  86-04.  Center 
for  Multie  .mate  Analesis.  February  1986. 

A  unified  approach  is  dee  eloped  for  measurement  and  analysis  of  diversity 
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in  data  with  a  factorial  structure,  i.e.,  classified  by  levels  of  a  number  of 
factors.  It  is  shown  that  the  total  diversity  in  given  data  can  be  decomposed 
into  several  components  each  representing  a  source  of  diversity,  by  using 
any  measure  of  diversity  satisfying  strong  convexity  properties.  It  is  shown 
that  the  usual  measures  of  diversity  like  Shannon's  entropy  have  only 
Limited  convexity  properties,  besides  being  not  applicable  to  quantitative 
data.  A  new  measure  of  diversity  called  quadratic  entropy  is  introduced 
and  its  applicability  to  both  qualitative  and  quantitative  data  is  discussed. 

68.  Taniguchi.  M..  Krishnaiah,  P.  R.,  and  Chao,  R.  Normalizing 
Transformations  of  Some  Statistics  of  Gaussian  ARMA  Processes. 
Technical  Report  No.  86-05,  Center  for  Multivariatee  Analvsis,  February 
1986. 

In  this  paper,  the  authors  investigate  Edgeworth  type  expansions  of  certain 
transformations  of  some  statistics  of  Gaussian  ARMA  processes.  They  also 
investigated  transformations  which  will  make  the  second  order  part  of  the 
Edgeworth  expansions  vanish.  Some  numerical  studies  are  made  and  they 
show  that  the  above  transformations  give  better  approximations  than  the 
usual  approximation. 

69.  Wu,  Yuehua,  On  Strongly  Consistent  Estimates  of  Regression  Coefficients 
When  the  Errors  are  not  Independently  and  Identically  Distributed, 
Technical  Report  #  86-06,  Center  for  Multivariate  Analysis,  March  1986. 

In  this  paper,  the  author  proposes  two  methods  of  estimation  of  the 
regression  coefficients  when  the  errors  are  not  disturbed  identically  and 
independently  and  are  of  nonzero  mean.  The  estimates  proved  in  this  paper 
are  shown  to  be  strongly  consistent  and  mean  square  consistent. 

70.  Chen,  X.R.  and  Zhao,  L.  C.,  Almost  sure  Lj-norm  Convergence  for  Data- 
based  Histogram  Density  Estimates,  Technical  Report  #  86-07,  Center  for 
Multivariate  Analysis,  March  1986. 

The  main  result  of  this  paper  is  summarized  in  Theorem  1.  which  states  that 
when  certain  conditions  of  a  general  nanire  are  satisfied,  the  data-based 
histogram  density  estimator  is  strongly  consistent  in  the  sense  that  the  mean 
absolute  deviation  of  the  estimator  and  the  density  function  converges  to 
zero  almost  surely  for  any  density  function,  as  the  sample  size  increase  to 
infinity. 

71.  Taniguchi.  M.  and  Krishnaiah.  P.R.,  Asymptotic  Distributions  of  Functions 
of  Eigenvalues  of  the  sample  Covariance  Matrix  and  Canonical  Correlation 
Matrix  in  Multivariate  Time  Senes",  Technical  Report  #  86-08,  Center  for 
Multivariate  Analysis.  March  1986. 

Let  X  1 1  ’.  .X  he  pxl  random  vectors  with  mean  zero. 

Put  S=-y;  .  When  X(t;,  it=l.  •  -  .n)  are  distributed  as  i.i.d.  NiO. 

X).  (i.e.  usuai  multivariate  analysis  ).  many  authors  have  investigated  the 
asymptotic  expansions  for  the  distributions  of  vanous  functions  of  the 
eigenvalues  of  S.  In  this  paper  we  will  extend  the  above  results  to  the  case 
when  (X_(tij  is  a  Gaussian  stationary  process.  Also  we  shall  dense  the 
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asymptotic  expansions  for  certain  functions  of  the  sample  canonical 
correlation  multivariate  time  senes.  Applications  of  some  of  the  results  in 
signal  processing  are  also  discussed. 

72.  Pukkila.  Tarmo  M.  and  Rao,  C.  Radhakrishna,  Pattern  Recognition  based 
on  scale  Invariant  Discriminant  Functions,  Technical  Report  #  86-09, 
Center  for  Multivariate  Analysis,  April  1986. 

Some  probability  models  for  classifying  individuals  as  belonging  to  one  of 
two  or  more  populations  using  scale  invariant  discriminant  functions  are 
considered.  The  investigation  is  motivated  by  practical  situations  where  the 
observed  data  on  an  individual  are  in  the  form  of  ratios  of  some  basic 
measurements  or  measurements  scaled  by  an  unknown  non-negative 
number.  The  probability  models  are  obtained  by  considering  a  p-vector 
random  variable  X  with  a  known  distribution  and  deriving  the  distribution 
of  the  random  vector  Y=[G(X)'l]X.  where  G(X)  is  a  non-negative  measure 
of  size  such  that  G(kX)=/.G(X)  for  /.>0.  Explicit  expressions  are  obtained 
for  the  densities  of  what  are  called  Angular  Gaussian,  Compositional 
Gaussian,  Type  1  and  Compositional  Gaussian,  type  2  distributions. 

73.  Krishnaiah,  P.R.,  On  Reduction  of  Dimensionality  Under  Multivariate 
Regression  Canonical  Correlation  Models,  Technical  Report  #  86-10, 
Center  for  Multivariate  Analysis,  April  1986. 

In  this  paper,  the  author  gives  a  review  of  the  literature  on  various 
techniques  for  determination  of  the  ranks  of  regression  matrix  and 
canonical  correlation  matrix.  Also,  methods  of  selection  of  important 
original  variables  under  multivariate  regression  canonical  ccTelation 
models  are  reviewed.  The  methods  reviewed  involve  not  only  tests  of 
hypothesis  but  also  model  selection  methods  based  ;pon  information 
theoretic  criteria. 

74.  Naik,  Dayanand  N..  Detection  of  Outliers  in  Multivariate  Linear  Regression 
Model,  Technical  Report  #  86-11,  Center  for  Multivariate  Analysis,  April 
1986. 

In  this  article  we  suggest  multivariate  kurtosis  measure  as  a  statistic  for 
detection  of  outliers  in  a  multivariate  linear  regression  model.  The  statistics 
has  some  local  optimal  properties. 

75.  Subramanyam.  K.  and  Rao,  M.  Bhaskara.  Extreme  Point  Methods  in  the 
Study  of  Classes  of  Bivariate  Distributions  and  Some  Applications  to 
Contingency  Tables.  Technical  Report  #  86-12.  Center  for  Multivariate 
Analysis.  April  1986. 

The  set  of  all  bivariate  positive  quadrant  distributions  is  neither  compact 
nor  convex.  But  the  set  of  all  bivariate  positive  quadrant  distributions  with 
fixed  marginals  is  a  convex  set.  These  convex  sets  tire  compact  in  the  case 
of  discrete  bivariate  distributions  if  the  marginals  have  finite  support.  A 
simple  method  to  enumerate  the  extreme  points  of  these  convex  sets  is 
given.  In  the  context  of  contingency  tables  for  testing  the  null  hypothesis 
independence  against  the  alternative  positive  quadrant  dependence  one  can 
use  the  method  of  extreme  point  analysis  to  compare  the  performance  of 
various  tests. 
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76.  Bai,  Z.D.  and  Fu.  J.C.,  Likelihood  Principal  and  Maximum  Likelihood 
Estimator  of  Location  Parameter  for  Cauchy  Distribution,  Technical  Report 
#  86-13,  Center  for  Multivariate  Analysis,  May  1986. 

In  the  literature  of  point  estimation,  Cauchy  distribution  with  location 
parameters  was  often  cited  as  an  example  for  the  failure  of  maximum 
likelihood  method  and  hence  the  failure  of  likelihood  principle  in  general. 
Contrary  to  the  above  notion,  we  proved,  even  in  this  case  that  the 
likelihood  equation  has  multiple  roots,  the  likelihood  estimator  (  the  global 
maximum)  remains  an  asymptotically  optimal  estimator  in  the  Bahadur 
sense. 

77.  Bai,  Z.D..Rao,  C.  Radhakrishna  and  Zhao,  L.C.,  Kernel  estimators  of 
density  function  of  directional  data.  Technical  Report  #  86-14,  Center  for 
Multivariate  Analysis,  May  1986. 

Let  X  be  a  unit  vector  random  variable  taking  values  on  a  k-dimensional 
sphere  Q  with  probability  density  function  f(x).  The  problem  considered  is 
one  of  estimating  f(x)  based  on  n  independent  observations  Xj,  ■  ■  •  .X^  on 
X.  The  proposed  estimator  is  of  the  form 

n 

fn(x)  =  (nhk-lrlC(h)J4K\(l-x'Xi)/h2],  xe  Q 

<=i 

where  K  is  a  kernel  function  defined  on  R+.  Conditions  are  imposed  on  K 
and  f  to  prove  pointwise  strong  consistency,  uniform  strong  consistency  and 
strong  Lpnorm  consistency  of  fn  as  an  estimator  off. 

78.  Pukkila,  Tarmo  M.  and  Krishnaiah.  P.R.,  On  the  use  of  autoregressive  order 
determination  criteria  in  univariate  white  noise  tests,  Technical  Report  # 
86-15,  Center  for  Multivariate  Analysis,  June  1986. 

In  this  paper  new  tests  of  white  noise  will  be  suggested..  All  of  the 
proposed  test  are  based  on  the  application  of  known  autoregressive  order 
determinating  criteria  such  as  AIC.BIC.  or  on  their  modifications.  Some  of 
the  proposed  tests  have,  among  other  things  the  property  that  the 
significance  level  of  a  test  approaches  to  zero  as  the  number  of  observations 
increases.  This  on  the  other  hand  means  that  if  this  kind  of  a  test  is  applied 
for  the  testing  whether  the  residual  series  of  a  fitted  model  is  white  noise, 
the  testing  leads  to  a  consistent  order  selection. 

79.  Pukkila,  Tarmo  M.  and  Krishnaiah.  P.R..  On  the  use  of  autoregressive  order 
determination  criteria  in  multivariate  white  noise  tests.  Technical  Report  # 
86-16.  Center  for  Multivariate  Analysis.  June  1986. 

In  this  paper  testing  the  hypothesis  of  multivariate  white  noise  is  seen  as  the 
selection  of  tnc  order  of  a  multivariate  autoregressive  model  for  the 
observed  time  senes.  Therefore  multivariate  white  noise  tests  can  be  carried 
out  by  apply  in;  autoregressive  order  determination  criteria  such  as  AlC. 
B1C.  etc.  I:  is  known  that  for  example  the  B1C  criterion  estimates 
consistently  the  order  of  an  autoregression.  An  order  determination 
criterion  with  this  property  leads  to  a  white  noise  test  with  the  significance 
level  approaching  to  zero  as  the  number  of  observerations  increases.  In  this 
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paper  the  order  of  autoregressive  moving  average  model  is  proposed  to  be 
determined  by  applying  this  kind  of  white  noise  test.  Then  the  resulting 
model  building  procedure  is  nothing  but  a  generalization  of  the  procedure 
proposed  by  Box  and  Jenkins. 

80.  Bai,  Z.D.  and  Yin,  Y.Q..  A  convergence  to  the  semicircle  law.  Technical 
Report  #  86-17,  Center  for  Multivariate  Analysis.  July  1986. 

In  this  article,  it  is  proved  that  the  spectral  distribution  of  the  random  matrix 

Ap=-±=(X/p-nIp) 

2 \np 

tends  to  the  semicircle  law  as  p— >  °°,  n  — >  ,o©  and  p/n— >0.  Here  X  = 
[X- ;  1  <i<p.  1  <j<n],  X-'s  are  iid  random  variables  with  EXj  j  =  0,  EXypl, 
EX|j<  oo. 

81.  Bai.  Z.D.,Krishnaiah,  P.R.  and  Yin.  Y.Q..  Inference  on  the 
Occurrence/Exposure  Rate  and  Simple  Rise  Rate,  Technical  Report  # 
86-18,  Center  for  Multivariate  Analysis,  August  1986. 

In  this  paper,  we  study  the  asymptotic  distributions  of  the  functions  of  the 
occurrence/exposure  rates  of  several  groups  of  patients  as  well  as  Berry- 
Esseen  bound  on  the  distribution  function  of  the  occurrence/exposure  rate. 
Asymptotic  distributions  of  functions  of  the  simple  risk  rates  are  also 
derived.  The  results  are  useful  in  not  only  medical  research  but  also  in  the 
area  of  reliability. 

82.  Das,  Rita  and  Sinha,  Bimal  K.,  Robust  optimum  invariant  tests  in  one-way 
unbalanced  and  tw'o-way  balanced  models.  Technical  Report  #  86-19. 
Center  for  Multivariate  Analysis.  August  1986. 

In  one-way  random  effects  unbalanced  model  the  locally  best  invariant  test 
for  the  quality  of  the  treatment  effects  is  derived.  Surprisingly,  this  is 
different  from  the  widelyused  familiar  F-test.  In  the  balanced  case, 
however,  the  tw'o  tests  coincide  and  represent  the  uniformly  most  powerful 
invariant  test.  For  two-way  random  effects  and  mixed  effects  balanced 
model,  the  uniformly  most  powerful  invariant  test  for  the  equality  of  the 
treatment  effects  is  derived  both  with  and  without  interaction,  and  shown  to 
be  equvalent  to  the  usual  F-test  under  fixed  effect  models.  The  optimum 
invariant  tests  derived  shown  not  to  depend  on  the  assumption  of  normality. 
Different  aspects  of  null,  nonnull  and  optimality  robustness  of  these  tests 
[Kariya  and  Sinha.  Annals  of  Statistics.  1985]  are  studied.  In  the 
unbalanced  two-way  models  however  unlike  in  the  fixed  efffects  model 
providing  a  l.'MP!  test,  both  random  and  mixed  effects  models  present  a 
difficulty  which  is  pointed  out. 

SC  Das.  Rita  and  Sinha.  Bimal  K..  Robust  optimum  invariant  tests  in 
covariance  structures  useful  in  linear  models.  Technical  Report  #  86-20. 
Center  for  Multivariate  Analysis.  August  1986. 

Necessary  and  sufficient  conditions  providing  structure  of  V  for  the  BLUE 
of  estimable  linear  parametric  functions  and  the  LRT  of  a  linear  testable 
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hypothesis  under  (yX0,  a2V)  are  well  known  in  the  literateure  (T.  Mathew 

and  P.  Bhimasankaram,  Sankhva  (A),  1983,221-225).  In  this  paper  we 
derive  robust  optimum  invariant  tests  of  such  structures  of  V  based  on  data 
generated  for  a  fixed  design  matrix  X.  Aspects  of  null,  nonnull  and 
optimality  robustness  of  the  proposed  tests  are  discussed. 

84.  Rao,  C.  Radhakrishna,  Applications  of  multivariate  analysis  m  signal 

detection:  some  recent  work.  Technical  Report  #  86-22,  Center  for 

Multivariate  Analysis.  August  1986. 

The  paper  contains  a  review  of  recent  research  in  the  detection  and 
estimation  of  signals.  Both  the  cases  of  fixed  and  random  signals  are 
considered.  The  effects  of  using  an  estimated  covanance  matrix  of  the  noise 
component  in  the  place  of  the  true  one  are  studied.  A  new  method  known 
as  the  efficient  detection  criterion  for  estimating  the  number  of  signals  is 
discussed. 

85.  Rao.  C.  Radhakrishna.  Some  recent  results  in  signal  detection.  Technical 
Report  #  86-23.  Center  for  Multivariate  Analysis,  September  1986. 

Some  recent  results  on  the  detection  and  estimtion  of  singnals  in  the 
presence  of  noise  are  discused.  An  exact  confidence  lower  bound  is 
obtained  for  the  discriminatory  power  of  an  estimated  linear  discriminant 
function  for  signal  detection.  Information  theoretic  criteria  are  suggested 
for  the  estimation  of  number  of  signals.  A  new  method  is  proposed  for 
determining  the  number  of  signals  and  estiomating  them  in  exponential 
signal  models. 

86.  Alzaid  A..  Lau.  K.S.,  Rao,  C.  Radhakrishna  and  Shanbhag.  D.N.,  Solution 
of  Deny  convolution  equation  restricted  to  half  line  via  random  approach. 
Technical  Report  #  86-24,  Center  for  Multivariate  Analysis.  September 
1986. 

A  general  solution  of  the  Deny  convolution  equation  restricted  to  half  line 
is  obtained  using  the  concepts  of  random  wralk  theory.  The  equation  in 
question  arises  in  several  places  in  applied  probability  such  as  in  queueing 
and  storage  theories  and  characterization  problems  of  probability 
distributions.  Some  of  the  important  applications  are  brief  lv  discussed. 

87.  Khatree.  Ravindra.  On  comparison  of  estimates  of  dispersion  using 
generalized  Pitman  nearness  criterion.  Technical  Report  #  86-25.  Center  for 
Multivariate  Analysis.  September  1986. 

In  this  paper,  estimates  of  dispersion  matrix  and  its  functions  are  compared 


based  on  generalized  Pitman  nearness  criterion.  Various  loss  functions  are 
considered  ;T-r  the  purpose.  Locally  superior  estimates  are  defined. 
Compari :  these  estimates  are  made  with  oilier  standard  ones.  It  is 
shown  that  w::h;:i  certain  classes,  defined  in  the  paper,  these  .ire  the  best 
estimate:  •  ,-:te:al;zed  Pitman  nearness  sense. 

88.  Bin.  /..  1)  \  mi  heller's  proof  on  a  law  of  iterated  logarithm.  Technical 
Report  ~  '■<  -!<  .  i  enter  ior  Multivariate  Analysis.  October  1986. 

In  194o.  hr.  er  popped  and  proved  a  famous  law  of  iterated  logarithm. 


But  unfortunately  his  proof  was  found  to  be  incorrect,  although  the  main 
result  was  true.  In  this  paper  the  author  gives  a  new  proof  of  the  main  result 
due  to  Feller. 

89.  Nishii.  R.  and  Krishnaiah.  P.R.,  On  the  moments  of  classical  estimates  of 
explanatory  variables  under  a  multivariate  calibration  model.  Technical 
Report  #  86-27.  Center  for  Multivariate  Analysis,  October  1986, 

In  the  area  of  regression  analysis,  the  problem  of  predicting  the  response 
variable  using  explanatory  variables  received  considerable  attention.  But 
situations  arise  when  it  is  of  interest  to  predict  the  explanatory  variable 
using  information  on  the  response  variables. 

90.  Yin,  Y.Q..  and  Krishnaiah,  P.R..  On  some  nonparametric  methods  for 
detection  on  the  number  of  signals.  Technical  Report  #  86-28,  Center  for 
Multivariate  Analysis,  October  1986. 

In  this  paper,  the  authors  considered  a  general  information  theoretic  criteria 
for  detection  of  the  number  of  signals  under  additive  model.  The  strong 
consistency  of  this  procedure  as  well  as  a  bound  on  the  probability  of 
wrong  detection  are  established.  Applications  of  this  method  to  detection  of 
the  number  of  signals  when  the  signals  are  coherent  is  also  discussed. 

91.  Bai.  Z.D.  .  Silverstein,  Jack  W.  and  Yin,  Y.Q.  A  note  on  the  largest 
eigenvalue  of  a  large  dimensional  sample  covariance  matrix.  Technical 
Report  #  86-29,  Center  for  Multivariate  Analysis,  Octooer  1986. 

Let  {Vjj;  i.j  =  1.2,  •  ■  •  }  be  a  family  of  i.i.d.  random  variables  with  E(vn4)= 
=».  For  positive  integers  p.n  with  p  =  p(n)  and  p/n— >  y  >  0  as  n  let 
Mn  =  0/n)VnVnT  where  Vn  =  (v^  <  and  let  X^fn)  denote  the  largest 

i  YSn 

eigenvalue  of  M^.  It  is  shown  that  lim  ^.max(n)  =°°  a.s.  This  result  verifies 

n 

the  boundedness  of  E(v  j  t4)  to  be  the  weakest  condition  known  to  assure  the 
almost  sure  convergence  of  Xmax(n)  for  a  class  of  sample  covariance 
matrices. 

92.  Kariya.  T.  and  Bimal  K.  Sinha.  and  Giri.  N.C..  Robustness  of  t-test. 
Technical  Report  #  86-31.  Center  for  Multivariate  Analysis.  October  1986. 

In  this  paper  we  establish  the  optimality  robustness  of  Student's  t-test 
mainly  without  invariance.  This  generalizes  some  well-known  results  of 
Lehmann  and  Stein  ( 1949).  and  Kariya  and  Eaton  1 1977). 

93.  Das.  Rita  and  Sinha.  Bimal  K..  Robust  optimum  invariant  tests  for  random 
MANOVA  models.  Technical  Report  #  86-32.  Center  for  Multivariate 
Analysis.  October  1986. 

Consider  tiie  canonical  form  MANOVA  setup  with  X:  nxp  =(.Y,'.AY..\Y)'  = 
(A/, h'  -  E.  X,  :  n.xp,  i=  1.2.3.  Mj  :  n,xp.  i=1.2.  nj  +  n,  -n,  =  n.  n,  > 
p.  where  E  i>  a  i.tndom  error  matrix  with  location  0  and  unknown  scale 
matrix  X  >0  ip.d.i.  Assume  .  unlike  in  the  usual  sense,  that  Mj  is  random 

with  location  0  and  scale  matrix  0|~X.  (}]“.  X  being  unknown.  For  testing 
Ho  :  G j ■  =  0  versus  H ]  :  a."  >  0  under  a  left  orthogonally  invariant 


distribution  of  X.  it  is  shown  that  when  either  n2  =  0  or  M-,  fixed  if  n2  >  0 
the  trace  test  of  Filial  0955)  is  UMPI  if  min  (nj,  p)  =  1  and  LBI  if  min  (n;. 
p)  >1.  The  test  is  null,  nonnull  and  optimality  robust  (Kariya  and 
Sinha(1985)).  However,  such  a  result  does  not  hold  if  n2  >  0  and  M2 
random. 

94.  Bai,  Z.D..Krishnaiah.  P.R.,  Rao,  C.  Radhakrishna,  Reddy,  P.S.,  Sun,  Y.N. 
and  Zhao.L.C.,  Reconstrucnon  of  the  left  ventricle  from  two  orthogonal 
projections.  Technical  Report  #  86-33.  Center  for  Multivariate  Analysis. 
November  1986. 

A  new  method  for  reconstruction  of  the  shape  of  the  left  ventricle  from  the 
biplane  angiocardiograms  is  proposed.  This  approach  utilizes  a  pair  of 
orthogonal  X-ray  projection  images.  The  shape  of  the  ventricle  is 
reconstructed  by  dividing  these  projection  images  into  parallel  slices  and 
then  being  processed  slice  by  slice  stepwise.  Each  corresponding  pair  of 
slices  from  two  one-dimensional  projection  profiles  which  are  used  to 
reconstruct  a  cross  section  of  the  ventricle. 

95.  Mathew.  Thomas  and  Sinha.  Bimal  Kumar.  Optimum  tests  for  fixed  effects 
and  variance  components  in  balanced  models.  Technical  Report  #  86-34, 
Center  for  Multivariate  Analysis.  November  1986. 

For  any  ANOVA  model  with  balanced  data  involving  both  fixed  and 
random  effects.  UMPU  and  UMPI  tests  are  derived  for  the  significance  of  a 
fixed  effect  or  a  variance  component,  under  the  assumption  of  normality  of 
random  effects.  The  tests  coincide  with  usual  F-tests.  Robustness  of  the 
UMPI  test  against  suitable  deviations  from  normality  is  established. 

96.  Sinha.  Bimal  K..  Inadmissibility  of  the  best  equivariant  estimators  of  the 
variance-covariance  matrix,  the  precision  matrix,  and  the  generalized 
variance:  A  survey.  Technical  Report  #  86-35.  Center  for  Multivariate 
Analysis.  November  1986. 

Based  on  a  data  matrix  X=(Xj,  ■  •  .Xk):  pxk  with  independent  columns 

X;  -  N  (Cj.  Z).  and  an  independent  p  x  p  Wishart  matrix  S  ~  Wp(n.I). 
procedures  to  obtain  estimators  dominating  the  best  equivariant  estimators 
of  Z.Z~!  and  El  under  various  loss  functions  are  reviewed. 

97.  Sinha.  Bimal  K.  and  Sen.  Pranab  Kumar.  Nonparametric  estimation  of  the 
generalized  variance.  Technical  Report  #  86-36.  Center  for  Multivariate 
Analysis.  November  1986. 

For  multivariate  distributions  with  finite  second  order  moments. a 
nonparametric  symmetric,  unbiased  estimator  of  the  generalized  variance  is 
considered  and  it  is  shown  to  be  (nonparametric)  optimal  for  the  class  of 
distribution  >  ha\  ing  finite  fourth  order  moments.  A  jackknife  version  of  the 
sample  generalized  variance  is  also  considered  as  a  contender;  it  is 
computationally  more  convenient  and  asymptotically  equivalent  to  the 
former.  It  i>  also  shown  that  the  second  estimator  performs  quite  well  (in 
large  sampiei  relative  to  the  optimal  normal  theory  estimators  under  several 
loss  functions. 
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98.  Bai,  Z.D.,  Krishnaiah.  P.R.  and  Zhao.  L.C..  On  simultaneous  estimation  of 
the  number  of  signals  and  frequencies  under  a  model  with  multiple 
sinusoids.  Technical  Report  #  86-37.  Center  for  Multivariate  Analysis, 
December  1986. 

In  this  paper,  the  authors  considered  the  problem  of  estimation  of  the 
frequencies  and  the  number  of  signals  under  a  signal  processing  model  with 
multiple  sinusoids.  The  frequencies  are  estimated  with  eigenvarianon  linear 
prediction  method.  The  number  of  signals  is  estimated  with  an  information 
theoretic  criterion.  The  strong  consistency  of  the  estimates  of  the 
frequencies  and  the  number  of  signals  is  also  established.  Also,  a 
modification  of  forward  backward  linear  prediction  method  is  suggested  to 
yield  consistent  estimators  of  the  frequencies. 

99.  Bai,  Z.D..  Krishnaiah.  P.R.  and  Zhao.  L.C..  On  rate  of  convergence  of 
equivariation  linear  prediction  estimates  of  the  number  of  signals  and 
frequencies  of  multiple  sinusoids.  Technical  Report  #  86-38.  Center  for 
Multivariate  Analysis.  December  1986. 

In  a  companion  paper,  the  authors  (see  Bai.  Krishnaiah  and  Zhao(1986)) 
considered  the  problem  of  estimation  of  the  number  of  signals  and  the 
frequencies  of  these  signals  under  a  signal  processing  model  with  multiple 
sinusoids.  The  number  of  signals  was  estimated  by  using  an  information 
theoretic  criterion.  They  have  also  established  the  strong  consistency  of 
their  estimates.  In  this  paper,  we  established  the  rates  of  convergence  of  the 
above  estimates  of  the  number  of  signals  and  frequencies. 

100.  Khatri.  C.G.,  Rao.  C.  Radhakrishna.Schaafsma.  W  .Steememan.  A.G.M. 
and  Van  Vark  G.N.  Inference  about  the  performance  of  Fisher's  linear 
discriminant  function  with  application  to  signal  detection.  Technical  Report 
#  86-39,  Center  for  Multivariate  Analysis.  December  1986. 

We  consider  the  problem  of  assessing  the  discriminatory  power  of  a  sample 
based  discriminant  function  for  classifying  a  p-vector  observation  as  arising 
from  one  or  the  other  of  two  alternative  distributions  (populations).  An 
important  application  of  tins  problem  arises  in  signal  detection  where  a 
decision  has  to  be  taken  whether  a  received  message  in  the  form  of  p-vector 
observation  is  pure  noise  or  signal  plus  noise.  When  the  signal  to  be 
detected  is  a  known  vector  5  and  the  noise  variance-covariance  matrix  is  a 
known  matrix  Z.  the  classification  is  done  through  Fisher's  linear 
discriminant  function  which  has  the  maximum  discriminatory  power.  If  5  or 
Z  or  both  are  unknown  we  estimate  the  LDF  by  substuting  for  6  and  Z  their 
estimates  based  .*u  previous  correctly  classified  samples.  However,  there  is 
loss  in  discriminatory  power  in  using  an  estimated  l.DF.  and  the  paper 
focuses  on  the  estimation  or  prediction  of  the  actual  discriminatory  power. 

The  disci ut. mat- -ry  power  ot  any  linear  discriminant  function  (LDF)  is 
measured  by  w :  >  Mtown  as  signal  to  noise  ratio,  which  is  the  difference 

m  expecie  !  ;.o  ot  the  l.DF  under  the  two  hypotheses  expressed  in 

standard  de\  units  In  order  to  estimate  signal  to  noise  ratio  we  use 

parametric  an.:  lease  one  out  methods. 

101.  Rao.  C.  Radhasrishna.  Prediction  of  future  observations  in  growth  curve 
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models.  Technical  Report  #  86-40,  Center  for  Multivariate  Analysis. 
December  1986. 

The  problem  of  predicting  a  future  measurement  on  an  individual  given  the 
past  measurements  is  discussed  under  non-parametic  and  parametric  growth 
models.  The  eficiencies  of  different  methods  of  prediction  are  assessed  by 
cross-validation  (CV)  or  leave-one-out  (LOO)  technique  in  each  of  three 
data  sets  and  the  results  are  compared.  Under  non-parametnc  models,  direct 
and  inverse  regression  methods  of  prediction  are  discribed  and  their  relative 
advantages  and  disadvantages  are  discussed.  Under  paramenc  models 
polynomial  and  factor  analytic  type  growth  curves  are  conssidered. 
Bayesian  and  empirical  Bayesian  methods  are  used  to  deal  with  unknown 
parameters. 

A  general  finding  is  that  much  of  the  information  for  forecasting  is 
contained  in  the  immediate  past  few  observations  or  a  few  summary 
statistics  based  on  past  data.  A  number  of  data  reduction  methods  are 
suggested  and  analyses  based  on  them  are  described.  The  usefulness  of  the 
LOO  technique  in  model  selection  is  demonstrated. 

102.  Tauxe.  W.N.  .Klein,  Hervert  A.  .Bagchi,  Ansuman  .Krishnaiah,  P.R. 
Kundu,  Debasis  and  Tepe,  Paul,  Clinical  evaluation  of  the  filtration 
fraction:  a  multivariate  statistical  analysis,  Technical  Report  #  86-41, 
Center  for  Multivariate  Analysis,  December  1986. 

Using  multivariate  statistical  technique,  we  have  separated  glomerular 
filtration  fraction  (FF)  data  into  various  clinical  groups.  In  general 
population  refered  for  renal  function  evaluation,  we  found  that  the  FF  did 
not  vary  with  gomerular  filtration  rate(GFR),  but  varied  inversely  with 
ERPF.  A  group  of  hypertensive  patients  was  identified  to  lie  sigmficandy 
above  the  FF/ERPF  regression  line  for  the  whole  group.  Another  group  of 
two-kidneyed  normal  subjects  was  identified  to  lie  below  this  regression 
line.  FF  of  a  group  of  one-kidneyed  normally  functioning  regression 
subjects  was  found  also  to  lie  similarly  disposed  below  this  regression  line 
but  to  have  significanlv  higher  FF  than  the  two-kidneyed  normals 
Coefficients  of  ellipses  that  defined  99^  confidence  regions  for  the  loci  of 
the  means  of  these  groups  were  derived. 

103.  Nishii.  R.  .  Bai.  Z.D.  .Krishnaiah.  P.R..  Strong  consistency  of  certain 
information  theoretic  criteria  for  model  selecnon  in  calibration, 
discriminant  analysis  and  canonical  correlation  analysis.  Technical  Report  # 
86-42  ,  Center  for  Multivariate  Analysis.  December  1986. 

In  this  paper,  the  authors  show  that  the  criteria  for  model  selection  based 
upon  efficient  detection  iBD)  criteria  are  consistent  for  certain  problems  in 
multivariate  calibration,  discriminant  analysis  and  canonical  correlation 
analysis.  I'hese  resuits  w  ill  be  proved  under  mild  conditions  on  the 
underlying  distribution. 

104.  Subraniunyam.  K  and  Rao.  Nl.  Bhaskara.  Analysis  of  total  positivity  of 
order  two  in  2xn  ordinal  contingency  tables.  Technical  Report  #  86-43  . 
Center  for  Multivariate  Analysis.  December  1986. 

The  set  of  nil  bivariate  probability  distributions  with  support  contained  in 
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{(i,j);  i=l,2  and  j=1.2.  ■  ,n{  which  are  totally  positive  of  order  two  is 

shown  to  be  a  convex  set  under  some  conditions  on  one  of  the  marginal 
distributions.  The  exteme  points  of  the  compact  convex  set  are  explicitly 
enumerated.  Using  the  structure  of  this  convex  set.  we  show  that  the  power 
function  of  any  test  for  testing  the  hypothesis  of  independence  against  the 
hypothesis  of  strict  total  positivity  of  order  two  in  2xn  ordinal  contingency 
tables  has  a  simple  form  in  terms  of  the  exteme  pooints.  A  numerical 
illustration  is  provided. 

105.  Shanbhag.  D.N.  and  Kotz.  S..  Some  new  approaches  to  multivariate 
probability  distributions.  Technical  Report  #  86-44.  Center  for  Multivariate 
Analysis.  December  1986. 

In  this  paper  we  extend  and  generalize  to  the  multivariate  set-up  our  earlier 
investigation  related  to  expected  remaining  life  funcnons  and  general 
hazard  measures  including  representations  and  stability  theorems  for 
arbitrary  probability  distributions  in  terms  of  these  concepts. 

106.  Niroomand  Chapeh.  H.A.  and  Rao,  M.  Bhaskara.  On  the  least  squares 
estimator  in  moving  average  models  of  order  one.  Technical  Report  # 
86-45.  Center  for  Multivariate  Analysis.  December  1986. 

A  simple  expression  is  derived  in  this  paper  for  the  error  sum  of  squares  in 
the  context  of  moving  average  models  of  order  one.  A  computer  program 
is  developed  to  estimate  the  parameter  of  a  moving  average  model  of  order 
one  based  on  the  method  of  least  squares. 

107.  Niroomand  Chapeh.  H.A.  ,Rao.  M  Bhaskara  .  A  new  method  of  estimation 
in  a  moving  average  model  of  order  one.  Technical  Report  #  86-46. 
December  1986. 

The  exact  likelihood  of  the  data  coming  from  a  moving  average  model  of 
order  one  is  complicated.  In  this  paper,  we  propose  a  method  of  estimation 
of  the  parameters  of  a  moving  average  model  of  order  one  based  on  the 
approximate  likelihood  of  the  data  and  on  the  simulation  of  a  pair  of 
random  variables.  Some  comparisons  were  made  of  this  with  some  weLl 
known  methods  for  moderate  sample  sizes.  A  computer  program  is 
appended  which  is  helpful  in  using  this  method. 

108.  Krishnaiah.  P  R.  .Miao.  B.Q.  and  Zhao.  L.C.,  On  detection  of  change 
points  using  mean  vector.  Technical  Report  ~  86-47  .  Center  for 
Multivariate  AnnKsis.  December  1986. 

In  this  paper,  the  authors  consider  the  problem  of  change  points  within  the 


framework  .  ■;  model  selection  procedures  using  information  theoretic 
criteria.  I'1'.  •  authors  proposed  procedures  for  estimation  of  the  locations  of 
ct'am  c  p.  ••  a:  d  the  number  of  change  points  The  strong  consistency  of 
these  j  .  :•  I  : -  ,t!-o  established.  Also,  the  problem  of  change  points  is 
disciis'cc  "  the  'ramework  ot  the  simultaneous  test  procedures. 

109.  Chen.  X  k  r. u.ih.  P.R.  and  Liang.  W.Q..  Estimation  of  multivariate 


btnar\  der.sitv  i.,-:ie  orthonormal  functions.  Technical  Report  #  86-48  . 
Center  ti  :  M-X.\ .mate  Analvsis.  December  1986. 


In  this  paper,  the  authors  studied  certain  properties  of  the  estimate  of  Liang 
and  Krishnaiah  ( J .  Multivariate  Anal.,  16.  162-172)  for  multivariate  binary 
density.  An  alternative  shrinkage  estimate  is  also  obtained.  The  above 
results  are  generalized  to  general  orthonormal  systems. 

1 10.  Chen.  X.R.  and  Krishnaiah,  P.R.  .  Test  of  linearity  in  general  regression 
models  .  Technical  Report  #  86-49.  Center  for  Multivariate  Analysis, 
December  1986. 

Let  (Xp  Y|).  i=l.  .n.  be  Lid.  samples  of  (X.  Y).  This  paper  proposes  a 

method  for  testing  the  linearity  of  the  regression  funcnon  E(Ylx  =  x).  The 
asymptotic  distribution  (under  null  hypothesis)  and  the  asymptotic  power  of 
the  test  are  determined.  Aiso.consistency  of  the  test  is  proved  under  mild 
conditions. 

111.  Chen.  X.R.  and  Krishnaiah.  P.R..  Estimation  and  testing  in  truncated  and 
non  truncated  linear  median-regression  models  .  Technical  Report  #  86-50  , 
Center  for  Multivariate  Analysis.  December  1986. 

Suppose  (Xj.V,).  i=l.  •  ,n  are  iid.  samples  of  (X.Y).  Instead  of  Yj.  we  can 

only  observe  Y,  =  maxiYj.O).  Denote  by  m(x)  the  median-regression 

function  of  Ywith  respect  to  X.  This  paper  discuss  the  estimation  of  m(x) 
when  it  is  assumed  that  /nU)=a+P'.t  for  some  a.(3.  The  consistency  and 
asymptotic  normality  of  the  estimators  (of  a  and  (3)  are  established.  Also,  a 
method  is  given  to  test  the  linearity  of  the  regression  function  m(x). 

1 12.  Nishn.  Rvuei,  Maximum  likelihood  principle  and  model  selection  when  the 
true  model  is  unspecified.  Technical  Report  #  87-01.  Center  for 
Multivariate  Analysis.  February  1987. 

Suppose  independent  observations  come  from  an  unspecified  distribution. 
Then  we  consider  the  maximum  likelihood  based  on  a  specified  parametric 
family  by  which  we  can  approximate  the  true  distribution  well.  We 
examine  the  asymptotic  properties  of  the  Quasi  likelihood  estimate  and  of 
the  quasi-maximum  likelihood.  These  results  wall  be  applied  to  model 
selection  problem. 

1  1?.  Nishii,  Ryuei.  Identities  and  inequalities  for  the  expected  values  of  random 
matrices  and  Loewner's  Theory.  Technical  Report  4  87-02.  Center  for 
Multivariate  Analysis.  March  1987. 

We  consider  the  identity  for  the  expected  values  of  the  quadratic  forms  of 
random  matrices.  As  a  special  case  the  identity  yields  the  matrix-valued 
inequality  oi  the  Schwarz  type.  The  most  important  application  is  for 
Loewner's  theorc  m  the  context  of  the  operator-monotone  function,  and 
mam  ir.equ.daies  h.i.ed  on  the  concavity  or  convexity  of  the  matrix-valued 
functions  .ev  dero.  cd  Further  generalization  of  positive  semi-definite 
matrices  u  f  1  'v  made. 

1  14.  Bai.  Z.  D  .  Krishnaiah.  P.R.  and  Zhao.  L.C..  Multivariate  components  of 
covariance  model  m  unbalanced  case.  Technical  Report  #  87-03.  Center  tor 
Multivariate  Altaic su.  March  1987. 
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In  this  paper  we  are  interested  in  the  determination  of  the  rank  of  the 
random  effect  covariance  matrix  in  an  unbalanced  components  of 
covariance  model,  and  obtain  some  strongly  consistent  estimates  of  the  rank 
based  on  eigenstructure  methods. 

115.  Chen.  X.,  Necessary  and  Sufficient  Conditions  for  the  Existence  of 
Consistent  Estimates  in  One-dimensional  Linear  Models  with  Truncated 
Error  Distribution,  Technical  Report  #87-04.  Center  for  Multivariate 
Analysis.  March  1987. 

Consider  the  linear  regression  model  with  one  independent  variable  x:  Y1  = 

a  +  (iXj  +  e  .  i  =  1 .  n.  where  a.  j3  are  unknown  parameters,  fx,}  is  a 

sequence  of  known  real  numbers,  and  |e()  is  a  sequence  of  nd.  random 
errors,  with  Ee,  =  0.  0  <  Var  e,  =  c~  <  a2  unknown.  Suppose  that  the 
probability  density  f(x)  of  et  vanishes  outside  a  bounded  interval 
and  on  [Cfj.C-i],  f(x)  is  positive  and  satisfies  the  Lipschitz  condition,  it  is 
shown  that  the  necessary  and  sufficient  condition  for  the  existence  of  a 
consistent  estimate  of  a.  (1  is  X^/H^  ->  0  and  1-^  ->  0  respectively,  where  ^ 

=  1/n^jL]  xl  and  The  case  in  which  medley  =  0  is  also 

considered,  resulting  in  the  same  necessary  and  sufficient  conditions. 

116.  Bai,  Z.D.  and  Yin.  Y.Q-  Necessary  and  Sufficient  Conditions  for  Almost 
Sure  Convergence  of  the  Largest  Eigenvalue  of  Wigner  Matrix.  Technical 
Report  #87-05.  Center  for  Multivariate  .Analysis,  Apnl  1987. 

Let  W  =  (X1J:12<i.j<°«)  be  an  infinite  matrix.  Suppose  W  is  symmetric, 
entries  on  the  diagonal  are  iid,  entries  off  the  diagonal  are  lid.  and  they  are 
independent.  Then  it  is  proved  that  the  necessary  and  sufficient  conditions 
for  >.max(l/\n  Wn)  -^a.  a.s..  are  (1)  E(X*])2  <  «.  (2)  EXj;  <  «>.  (?)  EX,s 
<  0.  (4)  a  =  2(7.  <72  =  EXt-,. 


17.  Bai.  Z.D.  .  Chen.  X.R.  .  Necessary  and  sufficient  conditions  for  the 
convergence  of  integrated  and  mean  integrated  P-th  order  error  of  the 
kernel  density  estimates.  Technical  Report  #  87-06.  Center  for  Multivariate 
Analysis.  Apnl  1987. 

Let  fjx)  =  (nh£)  'I^KUx-Xii/h,,)  be  the  kernel  estimate  of  a  d- 
dimentional  density  f.  where  the  kernel  function  K  is  a  probability  density 

on  Ra,  hn  >  0  is  a  constant,  and  Xj.  Xs . Xj.  are  iid.  samples  drawn  from 

the  population  with  density  f.  Let  p  2  1  and  define  m  =  |  !f,.(x)  - 
f{ x )IPdx .  Mrn  =  Etm  ).  It  is  shown  that  the  necessary  and  sufficient 
condition  for  any  one  of  the  following  two  assertions  lim^m,..,  =  0.  a.s  . 
lim  .  .  M.,„  =  0.  is  that  the  four  conditions:  |  f  "  \  'd.\  <  j  kl’txidx  < 
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1  18.  Chen.  X.R.  and  Zhao.  L.C.  .  Necessary  and  sufficient  conditions  for  the 
convergence  ot  integrated  and  me.ui  integrated  r-th  order  error  of  histogram 
density  estimates.  Technical  Report  #87-07.  Center  for  Multivariate 
Analysis.  .April  1087. 


Suppose  Xj.  . are  samples  drawn  from  a  population  with  density 
function  f.  and  fn(x)  =  fn(x:Xj,  -  .X^  is  an  estimate  of  f(x).  Denote  by 

m^  =  J  I  fn(x)  -  f(x)  lrdx  and  M^  =  Efm^)  the  integrated  r-th  Order  Error 
and  Mean  Integrated  r-th  Order  Error  of  fn  for  some  r>  1  (when  p=2.  they 
are  the  familiar  and  widely  studied  ISE  and  M3SE).  In  this  paper  the  same 
necessary  and  sufficient  condition  for  lim  =  0  and  lim  =  0  a.s.  is 

obtained  when  fn(x)  is  the  ordinary  histogram  estimator. 

119.  Bai.  Z.D. .Krishnaiah.  P.R.,  Rao.  C.R.  .  Sun.  Y.N.  and  Zhao.  L.C.. 
Reconstruction  of  the  shape  and  size  of  objects  from  two  orthogonal 
projections.  Technical  Report  #  87-08,  Center  for  Multivariate  Analysis. 
April  1987. 

Suppose  that  V  is  a  family  of  regular  regions  in  Rd  (d>2)  with  a  common 
constructive  profile  (CP),  i.e..  the  projection  along  xd-axis.  In  this  paper, 
the  authors  establish,  under  certain  restrictions  on  the  equal  divisor  (ED) 
the  one-to-one  correspondence  between  V  and  the  family  0  of  the  induced 
testing  profiles  i.e..  the  projection  along  Xj-axis  of  the  element  in  V.  This 

proves  that,  under  certain  restrictions,  a  regular  region  in  Rd  is  continuous 
function  of  its  profiles  (CP  and  TP),  which  theoretically  guarantees  the 
justification  of  the  reconstruction  of  the  left  ventrical  given  in  a  previous 
paper  of  ours. 

120.  Krishnaiah.  P.R.  and  Miao,  B.Q.  .  Control  chans  when  the  observations  are 
correlated.  Technical  Report  #  87-09,  Center  for  Multivanate  Analysis. 
May  1987. 

Traditionally,  control  chans  are  based  on  independently  normal  samples, 
but  in  practice  it  so  happens  that  the  samples  are  dependent.  In  this  review, 
dependent  samples  are  considered  as  ARMA  time  series.  Also, 
multidimensional  time  series  samples  are  discussed. 

121.  Rao.  C.  Radhakrishna  .  Estimation  in  linear  models  with  mixed  effects:  a 
unified  theory.  Technical  Report  #  87-10.  Center  for  Multivariate  Analysis. 
May  1987 

A  unified  theory  without  any  assumptions  on  the  ranks  of  the  matrices 
involved  is  developed  for  the  estimation  of  a  linear  function  of  unknown 
parameters  ((3)  and  hypothetical  variables  (c)  in  a  mixed  Gauss- Markoff 
linear  model.  Y=  X[i  +  L’c  +  t\  The  expressions  for  the  estimates  and  the 
mean  square  errors  of  the  estimates  depend  on  the  elements  of  .in  inverse 
partitioned  matrix  ax  in  the  case  of  the  Gausx-Markov  model  with  fixed 
effects  developed  earlier  bv  the  author  |Rao(  1971  )|. 

I  he  gene-.  1  the  i>  applied  to  several  special  prob'-mis  that  arise  in  social 
and  biolo  ■■  ..  x.  rercov  In  particular,  explicit  expressions  are  obtained  for 
ximuitait-.  mnu’i.'ii  or  prediction  in  several  similar  and  dissimilar 

linear  r  •  .i  without  concomitant  variables  and  m  special 

univariate  an.:  ariate  linear  models  in  genetic  selection. 

122.  Krishnaiah  1’  k  and  Miao.  B.Q.  .  Detecting  and  interval  estimation  about 
a  slope  change  pi  a  a  .  Technical  Report  #  87-11.  Center  for  Multivariate 
Analysis.  J une  1  os  ' 
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In  this  paper,  the  authors  consider  the  problem  of  change  points  using 
Gaussian  process.  The  distribution  of  the  statistic  to  estimate  a  change  pomt 
constructed  in  this  paper  can  be  approximated  by  the  first  type  of  extnmal 
distibution.  Based  on  this,  detection  and  interval  estimation  of  a  change 
point  in  various  situation  are  discussed. 


123.  Bai,  Z.D.  ,  Krishnaiah  P.R.  and  Zhao.  L.C.,  On  the  detection  of  arrival 
estimation.  Technical  Report  #  87-12.  Center  for  Multivariate  Analysis. 
June  1987. 


The  estimation  of  arrival  direction  is  important  task  in  signal  processing 
and  has  recently  received  considerable  attention  in  the  literature.  In  this 
paper,  the  authors  proposed  a  method  to  estimate  the  direction  of  arrival 
and  proved  the  strong  consistency  of  the  estimates  for  both  cases  in 
presence  of  white  noise  and  colored  noise. 


124.  Subramanvam  K.  and  Rao,  M.B.  ,  On  the  extreme  points  of  the  set  of  all 
2xn  bivariate  positive  quadrant  dependent  distributions  with  fixed 
marginals  and  some  applications  .  Technical  Report  #  87-13,  Center  for 
Multivariate  Analvsis.  June  1987. 


The  set  of  all  bivanate  distributions  with  support  contained  in  {(i.j);  i=l,2 
and  j=1.2.  •  •  •  ,n}  which  are  positive  quadrant  dependent  is  a  convex  set.  In 
this  paper,  an  algebraic  method  is  presented  for  the  enumeration  of  all 
exteme  points  of  this  convex  set.  Certain  measure  of  dependence,  including 
Kendall's  tau.  are  shown  to  be  affine  functions  on  this  convex  set.  This 
property  of  being  affine  helps  us  to  evaluate  the  asymptotic  power  of  tests 
based  on  these  measures  of  dependence  for  testing  the  hypothesis  of 
independence  against  strict  positive  quadrant  dependence. 


125.  Rao,  C.  Radhaknshna  .  Strategies  of  d^'a  analysis.  Technical  Report  # 
87-14.  Center  for  Multivariate  Analysis.  June  1987 


Styles  in  statistical  analysis  change  over  rime  while  the  object  of  "extracting 
all  the  information  from  data"  or  ’summarization  and  exposure"  remains 
the  same.  Statistics  has  not  yet  aged  into  a  stable  discipline  with  complete 
agreement  on  foundations.  Certain  methods  become  popular  at  one  time 
and  are  replaced  in  course  of  time  by  others  which  look  more  fashionable. 
In  spite  of  the  controversies,  the  statistical  methodology  and  fields  of 
applications  are  expanding.  The  computer  and  more  specifically  the 
availability  of  graphic  facilities  have  had  a  greater  impact  on  data  analysis. 
It  may  be  of  interest  to  briefly  review  some  historical  developments  in  data 
analvsis. 


126.  Wu.  Yuehu.i 
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Strong  consistency  of  estimation  of  number  of  regression 
ti  e  errors  are  independent  and  their  expectations  are  not 
cr  .  Technical  Report  -  X7-15.  (.'enter  for  Multivariate 
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random  errors  are  identically  distributed  and  have  zero  means,  smce  it  is 
sometimes  unrealistic.  As  a  compensation  for  this  relaxation,  we  assume 
the  errors  have  a  common  bounded  support  [aj.a2).  Under  certain 
conditions,  we  obtain  the  strongly  consistent  estimate  of  the  number  k  for 
which  (3  k  *  0  and  (3k+1  =  =  (3p  =  0.  by  using  the  information 

theoretical  criteria. 

127.  Bai,  Z.D.  Krishnaiah,  P.R.  and  Zhao.  L.C..  On  the  asymptotic  joint 
distributions  of  the  eigenvalues  of  random  matrices  which  arise  under 
components  of  covanance  model.  Technical  Report  #  87-16,  Center  for 
Mulrivanate  .Analysis.  June  1987. 

In  this  paper,  the  authors  derived  asymptotic  joint  distributions  of  the 
eigenvalues  of  some  random  matrices  which  arise  under  components  of 
covariance  model. 

128.  Bai.  Z.D.  .  Chen.  X.R.,  Krishnaiah.  P.R.  ,  Wu.  Y.H.  and  Zhao,  L.C. 
Strong  consistency  of  maximum  Likelihood  parameter  estimation  of 
superimposed  exponential  signals  in  noise.  Technical  Report  #  87-17, 
Center  for  Multivariate  Analysis,  June  1987. 

Consider  the  model  of  multiple  superimposed  exponential  signals  in 
additive  Gaussian  noise 

p 

Y,(t)  =  '^stJ}.l'  +  eJ(t),  r=0.1.  •  •  ,/t-l.  y=l.  ■  •  jV 

where  N  is  fixed  and  n— \  =  exp(\-l  C0j)).  co;6  [0,2rc),  i=l,  ■  ■  ■  ,p.  oij, 
s-  are  unknown  parameters  and  p  is  known.  Further,  e^t)  =  e^ft)  + 

\-l  ep(t),  and  e  j ( t).ej-»ft).  t=0. 1 .2,  •  •  •.  j=l,  ■  ■  ■  .N,  are  mumally 
independent  and  identically  distributed  real  random  variables  with  a 
common  distribution  N'(0,C“/2).  0  <  G~<  <*=.  o~  is  unknown.  It  is  shown  that 
if  0)j  *  CDj  when  i j  and  Z  j is-l  >  0  for  i=l.  .p,  then  the  Maximum 

Likelihood  Estimate  (/.j.  ■  •  .Ap)  is  strongly  consistent.  Moreover,  it  is 
shown  that  A.  converges  to  A,  with  ;ui  exponential  rate. 

12cf  Wu.  Y.H..  Strong  consistency  and  exponential  rate  of  the  minimum 
L.-norm  estimate  in  lineiu  regression  models  .  Technical  Report  #  87-18. 
Center  tor  Multivariate  Analysis.  June  1987. 

Let  y  =  (t.~\  '0  -  e..  i=l.  .n.  .be  a  linear  regression  mode!  where 
{ x, }  is  a  sequence  of  experimental  points,  i.e..  known  p-vectors.  (e  }  is  a 
sequence  ■  independent  random  errors,  with  medte.)=0.  i=  1.2. 

Define  the  :r.:n:i:u;m  1.,-norm  estimate  of  ia.jJ')'.  hv  (Ot  -.jj,,-)'.  to  be  chosen 
such  that 

>  \\  =  min  y  !\ u-.v/pi 

.  '  !  "  ,-l 

Under  quite  general  conditions  on  { x] }  and  te- }  the  strong  consistency  of 


the  minimum  Lj-norm  estimate  is  established.  Further,  under  an  additional 
condition  on  { } .  it  is  also  proved  that  for  any  given  e  >  0,  there  exist 
constant  C  >  0  not  depending  on  n.  such  that 

P{  II  &„-a  II2  +  II  (^-pll2  >  e2 !  <  exp{  -Cn } 

.  for  large  n. 

130.  Bai.  Z.D  Chen.  X.R.  Krishnaiah.  P  R.  and  Zhao,  L.C.  .  Asymptotic 
property  of  the  EVLP  estimation  for  superimposed  exponential  signals  in 
noise.  Technical  Report  #  87-19.  Center  for  Multivariate  Analysis.  Julv 
1987. 

This  paper  studies  the  model  of  superimposed  exponential  signals  in  noise: 

n 

Yj{t)  =  £  akf/.k‘  +  e,(f ),  r=0.1.  ■  ■  .n- 1,  j=  1.  ■  •  •  JV 

where  a. j.  •  •  ■  ./.  are  unknown  complex  parameters  with  module  1, 
/.  j,  are  unknown  complex  parameters  with  module  less  than  1, 

,/.p  are  assumed  distinct,  p  assumed  known  and  q  unknown. 
akj.k=l,--  ,p.  j=l.  ,N  are  unknown  complex  parameters.  ej(t), 

t=0. 1,  .n-1.  j=l.  ,N,  are  iid  complex  random  noise  variables  such 

that 

£e|(0).£le,(0)l2  =  b'.0<52<=c.£le1(0)l4<o° 

and  c2  is  unknown.  This  paper  gives: 

-  A  strong  consistent  estimate  of  q; 

-  Strong  consistent  estimates  of  .).q.62  and  lakjl,  k <  q; 

-  Limiting  distributions  for  some  of  these  estimates: 

-  A  proof  of  non-existence  of  consistent  estimates  for  /.k  and  akj,  k  >  q; 

-  A  discussion  of  the  case  that  N  — >  °° 

131.  Burbea  Jacob  and  Oiler.  Jose  M.  .  The  information  metric  for  univariate 
linear  elliptic  models  .  Technical  Report  #  87-20.  Center  for  Multivariate 
Analysis.  June  1987 

The  information  metric  associated  with  a  univariate  linear  elliptic  family  is 
shown  to  be.  Essentially,  the  Poincare’  hyperbolic  metric  on  a  half-space 
whose  geodesic  Rao  distance  is  an  increasing  hyperbolic  function  of  a 
modified  Mahalanohis  distance.  This  result  enables  us  to  construct  new 
statistical  tests  and  to  recover  earlier  results  as  special  cases. 

132.  Subramam.ini  K.  R..o.  M  Bhask.ira  .  Analysis  ot  regression  dependence  in 

2\n  hr. a: .::,>n>  and  some  applications  in  contingency  tables  . 
Technical  -  '"-21.  Center  for  Multivanate  Analysis.  July  1987 

Tlte  not:  ■  ■  p. e  regression  lependence  among  discrete  bivariate 

distnbutio:.'  .  -  examined  vis-a-vis  other  well-known  notions  of 
dependence  Eve  ::e:ne  points  of  the  set  of  till  bixanate  2xn  distributions 
m  which  Y  . '  :  -gression  dependent  on  X  are  enumerated.  Some 
applications  t.  n.c  .•}  b.y po theses  are  given. 


133.  Krishnaiah  P.R.  Miao,  B.Q.  and  Zhao,  L.C..  Local  likelihood  method  in 
the  problems  related  to  change  points  ,  Technical  Report  #  87-22,  Center 
for  Multivariate  Analysis.  June  1987. 

In  this  paper,  the  so-called  local  likelihood  method  is  suggested  for  solving 
the  change  point  problems  when  the  data  are  distributed  as  multivariate 
normal.  The  detection  procedures  proposed  not  only  provide  strongly 
consistent  estimates  for  the  number  and  location  of  the  change  points,  but 
also  simplify  significantly  the  computation. 

134.  Bai.  Z.D.  ,  Krishnaiah.  P.R.  and  Zhao.  L.C..  Variable  selection  in  logistic 
regression.  Technical  Report  #  87-23.  Center  for  Multivariate  Analvsis. 
July  1987. 

In  many  situations,  we  are  interested  in  selection  of  important  variables 
which  are  adequate  for  prediction  under  a  logistic  regression  model.  In  this 
paper,  some  selection  procedures  based  on  the  information  theoretic  criteria 
are  proposed,  and  these  procedures  are  proved  to  be  strongly  consistent. 

135.  Chen,  X.R  and  Yuehua,  Wu  .  Strong  consistency  of  M-estimates  for  the 
linear  model  ,  Technical  Report  #  87-24.  Center  for  Multivariate  Analvsis, 
July  1987. 

This  article  defines  the  M-estimate  for  the  linear  model  directly  from  the 
minimization  problem 

n 

Xp(y,-cx-P'X,)  =  min. 

i=i 

Suppose  that  (X^Y,).  ■  ,(Xn.Yrl),  are  i.i.d.  observations  of  a  random 

vector  (X.Y).  where  Y  is  one-dimensional  and  X  may  be  multidimensional. 
It  is  shown  that  the  M-estimates  (3n  defined  in  this  manner  converge  with 
probability  one  to  o^.Pq  respectively  ((cxq.Pq’)  is  the  true  parameter)  as 
n— under  verv  general  conditions  on  the  function  p  and  the  distribution 
of  (X.Y). 

136.  Chen.  Xiru  .  Tesnng  and  interval  estimation  in  a  change-point  model 
allowing  at  most  one  change  .  Technical  Report  ~  87-25.  Center  for 
Multivariate  .Analysis,  July  1987. 

This  paper  considers  the  simplest  model  of  change-point  in  which  at  most 
one  change  in  the  mean  may  occur.  Results  include: 

-  Introduction  of  a  test  for  the  null  hypothesis  that  "no  change  in  the 
mean  occurs",  and  the  limit  distribution  of  the  test-statistic. 

-  Approximate  calculation  of  the  power  of  the  test. 

-  Interval  estimation  of  the  position  of  change. 

-  Point  estimation  of  the  ]ump  at  the  point  of  ciiange  and  its  asymptotic 
distribution. 

-  Evaluation  of  the  bias  of  the  MI.E  of  error  variance. 

137.  Rao.  C.  Radhakrishna.  and  Shanbhag.  D.N..  Further  Extensions  of  the 
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Choquet-Deny  and  Deny  Theorems  With  Applications  In  Characterization 
Theory,  Technical  Report  #87-26,  Center  for  Multivariate  Analysis.  August 
1987." 

Davies  and  Shanbhag  (1987)  identified,  under  some  mild  conditions,  the 
solution  to  the  convolution  equation  H(x)  =  j  sH(x+y)p.(dy).  x  e  S  where  S 
is  a  Polish  Abelian  semigroup  with  zero  element,  H:  S— »  R+.  a 
nonnegative  continuous  function  and  pi.  is  a  measure.  A  variant  of  the  result 
in  the  case  where  H  is  bounded  and  p.  is  a  certain  bounded  signed  measure 
is  obtained.  This  provides  a  generalized  version  of  the  Choquet-Deny 
theorem  where  p.  is  considered  to  be  a  probability  measure  and  S  to  be  a 
group.  The  stability  of  the  solution  to  the  above  equation  is  examined  by 
adding  an  error  term.  Further,  solutions  to  simultaneous  equations  of  the 
above  type  are  considered.  The  results  of  this  paper  generalize  those 
obtained  by  Gu  and  Lau  (1984).  Lau  and  Rao  (1982),  Ramachandran  et  al 
(1987),  Shimuzu  (.1978)  and  others  for  S  =  R  or  R+. 

138.  Yin,  Y.Q.  Detection  of  the  Number,  Locations  and  Magnitudes  of  Jumps, 
Technical  Report  #87-27,  Center  for  Multivariate  Analysis.  August  1987. 

Consider  a  signal  x(t)  =  fit)  +  w(t),  0  Ite  t  <  1.  Here  the  noise  w(t)  is  an 
independent  process,  and  fit)  is  a  function  with  only  finitely  many  jumps, 
satisfies  a  Lipschitz'  condition  between  any  two  consecutive  jumps.  This 
paper  gives  an  algorithm  to  determine  the  number,  locations  and 
magnitudes  of  the  jumps  of  fit).  The  consistency  and  speeds  of 
convergence  are  obtained. 

139.  Yin.  Y.Q.  .An  Effective  Method  to  Estimate  the  Number  and  Locations  of 
Jumps.  Technical  Report  #87-28.  Center  for  Multivariate  Analysis.  August 
1987. 

Let  x(t)  =  fit)  +  w(t)  be  a  signal,  where  the  noise  w(t)  is  an  independent 
process.  Suppose  fit)  has  only  finitely  many  jumps,  satisfies  a  Lipschitz' 
condition  between  each  two  consecutive  jumps.  This  paper  gives  a  very 
simple  algorithm  to  determine  the  number  and  locations  of  the  jumps.  The 
consistency  is  established,  the  speed  of  convergence  is  obtained. 

140.  Zhao.  L.C.  and  Chen.  X.R.  Normal  Approximation  For  Finite-population 
L-statistics.  Technical  Report  #87-29.  Center  for  Multivariate  Analysis, 
August  1987. 

In  this  article,  the  general  central  limit  theorem  and  Berry-Esseen  bounds 
for  finite-population  l /-statistics  with  degree  m  are  established  under  very 
weak  conditions.  The  results  achieve  a  substantial  improvement  of  those  in 
Zhao  and  ( 'hen  i  !  487  >. 

141.  Zhao.  !.(’  Eii'hr.aiah.  P.R.  and  Chen.  X.R..  Almost  Sure  l  -norm 
Converge:!.. •  l  or  Data  based  Histogram  Density  Estimates.  Technical 
Report  -s  --1  Center  tor  Multivariate  Analysis.  .August  1487, 

Let  X  \  v  ■  ■  .!  samples  drawn  from  a  d-dimensional  distribution  with 
density  t.  Partition  the  space  into  a  union  of  disjoint  intervals  f  I,  = 
I(/.X | . X,  w  i;|;  the  tonn 


•  •  ■  Jtf1*) — °°<all<xi‘,)<b ^<00,/=  1 .  -  -  .d) 

Define  the  data-based  histogram  estimate  of  f(x)  based  on  this  partition  as 

fn(x)  =  The  number  of  X, . falling  into  /  n  times  the  volume  of  I/,  for 

x  e  //•  1  =  !-2. For  given  constant  r  >  1  we  obtain  the  sufficient  condition 
for  Lim^Jf^xl-flxlFdx  =  0.  The  results  give  substantial  improvements 
upon  the  existing  results. 

142.  Shen.  W.H.  and  Sinha,  B.K.  Adimssible  Bayes  Tests  For  Structural 
Relationship,  Technical  Report  #87-31.  Center  for  Multivariate  Analvsis. 
September  1987. 

It  is  an  open  problem  to  construct  a  test  for  structural  relationship  among 
the  mean  vectors  of  several  multivariate  normal  populations  with  know-n 
but  unequal  covariance  matrices.  In  this  paper,  a  class  of  admissible  Bayes 
tests  for  the  above  problem  is  derived.  As  a  byproduct,  in  the  special  case 
of  known  and  equal  covariance  matrices,  the  likelihood  ratio  test  of 
Rao(1973)  is  shown  to  be  admissible  Bayes. 

143.  Subramanyam.  K.  Some  Comments  On  Positive  Quadrant  Dependence  In 
Higher  Dimensions.  Technical  Report  #87-32.  Center  for  Multivariate 
Analysis.  September  1987. 

An  extreme  point  analysis  has  been  performed  on  two  natural  definitions  of 
positive  quadrant  dependence  of  three  random  variable.  This  analysis  helps 
us  to  understand  how  much  these  two  notions  of  dependence  are  different. 

144.  Krishnaiah.  P.R..  Miao.  B.Q.  and  Wong.  H.  Multidimensional  Control 
Chart  For  Correlated  Data.  Technical  Report  #87-33.  Center  for 
Multivariate  Analysis.  June  1987. 

rhe  authors  discuss  the  problem  of  multivariate  control  charts  by 
considering  observations  as  time  senes. 

145.  Ghosh.  M.  and  Sinha.  B.K.  Empirical  and  Hierarchical  Bayes  Competitors 
of  Preliminary  Test  Estimators  In  Two  Sample  Problems.  Technical  Report 
#87-34,  Center  for  Multivariate  Analysis.  September  1987. 

Suppose  in  a  laboratory,  say  Laboratory  I.  a  certain  instrument  is  designed 
to  measure  several  characteristics,  and  a  number  of  vector-valued 
meas'iremcnts  ls  recorded.  Our  objectives  to  estimate  the  unknown 
population  mean.  It  is  known,  however,  that  a  similar  instrument  is  used  in 
another  laboratory,  say  Laboratory  II  for  the  same  purpose,  and  a  number  of 
observations  is  recorded  trom  the  second  instrument  It  is  also  suspected  that 
the  two  population  means  are  equal,  in  which  case,  obsrvauons  recorded  in 
Laboi. time  II  possibly  be  used  eflectivelv  together  wuh  those  in 

Laboratory  I  mr  estimating  the  population  mean  of  the  first  instrument. 
Thus,  the  qmw-a  that  naturally  arises  is  whether  one  should  use  the 
sample  meat:  tto.u  Laboratory  I  or  the  pooled  mean  from  the  two 
laboratories 


In  problem 


to  type  what  is  normally  sought  is  a  compromise  estimator 
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which  leans  more  towards  the  pooled  sample  mean  when  the  null 
hypothesis  of  the  equality  of  the  two  population  mean  is  accepted,  and 
towards  the  sample  mean  from  Laboratory  I  which  when  such  a  hypothesis 
is  rejected. 

In  this  paper,  we  propose  instead  an  empirical  Baves(EB)  estimator  which 
achieves  the  intended  compromise.  Such  an  EB  estimator  is  quite  often  a 
weighted  average  of  the  pooled  mean  and  the  first  sample  mean.  Unlike  a 
subjective  Bayes  estimator,  the  EB  estimator  is  quite  robust  (with  respect  to 
its  frequentist  or  Bayesian  risk)  against  a  wide  class  of  priors. 

146.  Bai.  Z.D..  Chen,  X.R..  Wu,  Y.  and  Zhao.  L.C.  Asymptotic  Normality  of 
Minimum  Lj-norm  Estimates  In  Linear  Models,  Technical  Report  #87-35, 
Center  for  Multivariate  Analysis,  September  1987. 

Consider  the  standard  linear  model  Yj  =t'0  +  ej,  i  =  1 . n,...  where  x>,  x2, 

...  are  assumed  to  be  known  p-  vectors.  0  the  unknown  p-vector  of 
regression  coefficients,  and  ej,  e-,,  ...  the  independent  random  error 
sequence  each  having  a  median  zero.  Define  the  Minimum  Lj-Norm 

estimator  as  the  solution  of  the  minimization  problem  £{Li  I  Y;-xJ^  I  = 

inff£"_ll  Yj  -  xj0  I:  0  e  RH.  It  is  proved  in  this  paper  that  0^  is 
asymptotically  normal  under  very  weak  conditions.  In  particular,  the 
condition  imposed  on  { xj }  is  exactly  the  same  which  ensures  the 
asymptotic  normality  of  Least  Squares  estimate:  lLmn->~maxi<i<n 

xi(Xjn=lSj  ^  =  0. 

147.  Bhandary  M.  Robust  M-estimation  of  Dispersion  Matrix  With  a  Structure. 
Technical  Report  #87-36.  Center  for  Multivariate  Analvsis.  September 
1987. 

Marona  (1976)  considered  the  problem  of  robust  M-estimation  of  the 
dispersion  matrix  X  which  is  arbitrary.  In  this  paper  a  procedure  for  robust 
M-  estimation  of  the  dispersion  matrix  of  the  form  T  +  (7‘Ip.  where  T  is 

n.n.d.  of  rank  q(<p)  and  o~  >  0  are  unknowns,  is  given.  The  structure  of  the 
above  covariance  matrix  appears  in  the  area  of  signal  processing.  The 
procedure  is  an  algorithm  which  involves  iterative  solution  of  the 
covariance  matrix.  This  procedure  takes  care  of  outliers  in  the  data  and  in 
fact  reduces  the  effect  of  outliers.  Simulation  result  shows  that  the  iteration 
converges  after  some  steps.  The  consistency  and  asymptotic  normality  of 
the  estimator  follows  from  Huber's  ( 1967)  result. 

148.  Bhandary  M.  Detection  of  Outliers  in  Signal  Processing  and  Detection  of 
Number  of  Signals  in  Presence  of  Outliers.  Technical  Report  #87-37. 
Center  for  Multivariate  Analysis.  September  1987. 

Likelihood  ratio  tests  for  mean-slippage  outliersisi  and  Roy's  1 1953)  union- 
intersection  p-ixiple  for  dispersion-slippage  outliensi  are  applied  in  the 
signal  pro.  a-om:  data  Procedures  to  get  approximate  critical  values  of  the 
tests  are  go.  am  1  hen  likelihood  ratio  tests  and  intuitive  tests  are  suggested 
for  mean-si. ppage  and  dispersion-slippage  outliers  respectively  to  test  the 
rank  of  I  tor  the  covariance  matrix  iT  +  a-!)  in  the  presence  of  outliers. 


Finally  the  estimate  of  rank  of  T  is  obtained  through  information  cntenon 
developed  by  Zhao.  Krishnaiah  and  Bai  (1986a.b)  and  through  modified 
information  cntenon  proposed  by  the  author. 

149.  Bhandary  M.  Detection  of  the  Number  of  Signals  in  the  Presence  of  White- 
noise  in  Elliptically  Symmetric  Distnbution.  Technical  Report  #87-38. 
Center  for  Multivanate  Analysis.  September  1987. 

A  new  method  for  the  detection  of  the  number  of  signals  in  the  white-noise 
case  when  the  underlying  distnbution  is  elhpncaily  symmetric  is  proposed. 
This  new  method  suggests  an  algorithm  for  the  maximum  likelihood 
estimation  of  parameters  which  involves  iteranve  procedure.  The  estimate 
of  parameters  reduce  the  effect  of  outliers.  Finally  the  estimate  of  number 
of  signals  is  obtained  tiuough  modified  information  cntenon  proposed  by 
the  author  and  the  estimator  is  strongly  consistent. 

150.  Bhandary  M.  Detection  of  the  Number  of  Defecnve  Components  of  Outlier 
Observations  m  Signal  Processing  Data.  Technical  Report  #87-39.  Center 
for  Multivanate  Analysis.  September  1987. 

Methods  of  detection  of  number  of  defective  components  in  the  outlier 
observations  in  signal  processing  data  have  been  suggested.  The  problem  is 
solved  for  both  mean-slippage  and  dispersion-slippage  outliers  by  Roy's 
(1953)  and  by  modified  information  cntenon  proposed  by  the  author.  The 
esnmators  are  strongly  consistent. 

15 1.  Tamguchi.  M..  Zhao,  L.C..  Knshnaiah.  P.R.  and  Bai.  Z.D.  Statisncal 
Analysis  of  Dyadic  Stationary  Processes.  Technical  Repon  #87-40.  Center 
for  Mulnvanate  Analysis.  September  1987. 

In  this  paper  we  consider  a  multiple  dyadic  stationary  process  with  the 
Walsh  spectral  density  matrix  fQ t  a.).  where  (j)is  an  unknown  parameter 

vector.  We  define  a  quasi-maximum  likelihood  estimator  £)of  Q  and  give 
the  asymptotic  distnbution  of  funder  appropnate  conditions.  Then  we 
propose  an  information  criterion  which  determines  the  order  of  the  model, 
and  show  that  this  cntenon  gives  the  consistent  order  estimate.  As  for  a 
finite  order  dyadic  autoregressive  model,  we  propose  a  simpler  order 
determination  criterion,  and  discuss  its  asymptotic  properties  in  detail  . 
This  criterion  gives  strong  consistent  order  estimate.  Also  detections  of 
signals  for  dyadic  stationary  processes  will  be  discussed.  In  Section  P  we 
discuss  testing  whether  an  unknown  parameter  O  satisfies  a  linear 
restriction.  Then  we  give  the  asymptotic  distribution  of  the  likelihood  ratio 
criterion  under  the  null  hy  pothesis. 

152.  Bai.  Z  D..  Kr:si,:uiah.  P.R.  and  Zhao.  L.C.  On  Multiplicities  of  the 
Eigenvalues  Cot;  p. menus  of  Covariance  With  Applications  In  Signal 
Proce.s'im:.  I  echm.al  Report  -87-41.  Center  for  Multivariate  Analysis. 
Septetnbei  !'•' 

In  this  pa  pet.  authors  investigated  procedures  for  determm  ion  of  the 
multiplicities  the  eigenvalues  of  the  covariance  matrix  of  the  random 
effects  vect.-r  under  multivariate  random  effects  model. 
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1.  Knshnaiah,  P.  R.  "Multivariate  Gamma  Distribution"  Encyclopedia  of 
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2.  Knshnaiah  P.  R.  "Multivariate  Multiple  Comparisons,"  Encyclopedia  of 
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4.  Bar.  Z.  D.  "The  Strong  Consistency  of  Error  Probability  Estimates  in  NN 
Discrimination."  Chinese  Annals  of  Math,  Vol.  6.  No.  3.  299-308.  1985. 
(Technical  Report  #84-43) 
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Applications.  Vol.  7().  105-1 14,  1985.  (Technical  Report  #84-27) 

6.  Rao.  C.  Radhakrishna,  "The  Inefficiency  of  Least  Squares:  Extensions  of 
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Equations."  Probability  Theory  and  Related  Fields  Vol.  73,  539-553.  1986. 
(Technical  Report  #84-41) 
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